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Abstract 

We present new results on the Gross-Neveu model at finite temperature and at next-to-leading 
order in the expansion. In particular, a new expression is obtained for the effective potential 
which is explicitly invariant under renormalization group transformations. The model is used as 
a playground to investigate various features of field theory at finite temperature. For example 
we verify that, as expected from general arguments, the cancellation of ultraviolet divergences 
takes place at finite temperature without the need for introducing counterterms beyond those of 
zero-temperature. As well known, the discrete chiral symmetry of the 1+1 dimensional model is 
spontaneously broken at zero temperature and restored, in leading order, at some temperature T^, 
we find that the 1/A^ approximation breaks down for temperatures below Tc- As the temperature 
increases, the fiuctuations become eventually too large to be treated as corrections, and a Landau 
pole invalidates the calculation of the effective potential in the vicinity of its minimum. Beyond Tc, 
the 1/A^ expansion becomes again regular: it predicts that in leading order the system behaves as a 
free gas of massless fermions and that, at the next-to-leading order, it remains weakly interacting. 
In the limit of large temperature, the pressure coincides with that given by perturbation theory 
with a coupling constant defined at a scale of the order of the temperature, as expected from 
asymptotic freedom. 
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I. INTRODUCTION 



This paper presents a detailed and pedagogical study of the Gross-Neveu model [|T| at 
finite temperature at next-to-leading order in the expansion. As well known, the model 
describes a system of interacting fermions in one spatial dimension. It is renormalizable 
and asymptotically free, and, in the vacuum, exhibits chiral symmetry breaking. These 
properties, chiral symmetry breaking and asymptotic freedom, mimic those of Quantum 
Chromodynamics (QCD) and the Gross-Neveu model constitutes an ideal playground to 
study these questions in a much simpler context than that of non abelian gauge theories in 
four dimensions. Because of this connection, we shall often use the language of QCD in this 
paper and refer to the fermion as a "quark" , and to the fermion species as "flavors" . 

Aside from being perturbatively renormalizable, the model is also renormalizable in the 
1/A^ expansion for any dimension smaller than four P, ^, ^, |, ^. It is one of the main 
motivations of the present work to study the details of the renormalization in a scheme 
which is not restricted to perturbation theory. In particular, it is expected on general 
grounds 0, H, |^ that the short distance singularities leading to ultraviolet divergences 
are not affected by the temperature. That is, inflnities which occur in flnite temperature 
calculations can all be removed by the zero temperature counterterms. It is interesting to 
see explicitly, on a non trivial example, how these cancellations of inflnities take place. 

Because the calculations are simpler in one dimension we restrict ourselves here to this 
situation. But, our main interest is not the one- dimensional physics and we leave aside 
aspects of the model which are speciflc to one dimension, at the cost of obscuring perhaps 
the physical interpretation of some of our results. In particular, the 1/A^ expansion is built on 
fleld conflgurations which make the action an extremum, and in this paper we consider only 
static, uniform, such conflgurations. However the action is extremal also for conflgurations 
which are not uniform, namely conflgurations commonly called "kinks", which interpolate 
between the two degenerate minima of the leading order effective potential. The role of these 
kinks at flnite temperature has been qualitatively discussed for the Gross-Neveu model in 
Ref . |]10|] . The kinks have energies typically of order e ~ NMf where Mf is the fermion mass 
(arising from symmetry breaking), and their number is ~ Le~^/-^, where L is the length of 
the system. Their role depends then of the order of the two limits A^ oo, L —>■ oo. If one 
takes the limit N oo flrst, then the kinks play no role since their number is exponentially 
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small: this is the situation described by the mean field approximation, or the leading order 
in the expansion, and where symmetry breaking occurs. When going to the next- 
to- leading order in the expansion, one has to take the thermodynamic limit L ^ oo 
at fixed A^; in this case the entropy associated with the positions of the kinks eventually 
overcome the cost in energy for producing them. These configurations are then expected to 
dominate, preventing symmetry breaking at any non-zero temperature, in agreement with 
general results about infinite systems in one-dimension |Tl|. Note also that kinks play a 
major role in the exact solution of the model found recently [0. 

In our calculations, kinks are ignored. Then, in the leading order of the expansion, 
chiral symmetry is restored at some finite temperature T^. But the calculation of the cor- 
rections of order reveals that the expansion breaks down as one approaches T^: The 
fiuctuations become eventually too large to be treated as corrections, which may be related 
to the existence of a Landau pole invalidating the calculation of the effective potential for 
small values of the field. It is unclear to us whether such difficulties with the expansion 
are consequences of the one dimensional character of the system, and in particular whether 
they could be cured for instance by taking kinks into account. Since treating kinks explicitly 
would represent a major effort beyond the scope of this paper, we leave this question open 
and focus on another interesting regime, that of high temperature. In this regime, because of 
asymptotic freedom, one expects the system to become a weakly interacting gas of fermions, 
somewhat similar to the quark-gluon plasma of QCD. As we shall see, verifying explicitly 
how such properties emerge in the expansion is instructive. 

There has been many studies of the Gross-Neveu model, or of the related Nambu Jona- 
Lasinio model |jl3|, at finite temperature. Most of these studies however concentrated on 



the mean field physics and chiral symmetry breaking [0, |T5|, |T6|, [T^, 0, |T9|, |20|, More 
recently, corrections where calculated in a version of the Gross-Neveu model with con- 
tinuous symmetry, focusing on the dominant role of the soft pion modes at low temperature 



2% , p3| . Investigations of the effect of the corrections in the three-dimensional Nambu 



Jona-Lasionio model were also presented in Refs. However, none of these studies 

provide systematic answers to the set of questions that we address in the present work. 

The outline of this paper is as follows. The next section is a general introduction to the 
model; we recall the construction of the effective potential at finite temperature and of its 
expansion. In Sect. |ITJ we calculate explicitly the zero temperature effective potential 
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at next-to- leading order in the expansion and carry out completely its renormaliza- 
tion. Although there exists many calculations of the effective potential at zero temperature 
p^ , p8| , p9[| , the expressions that we obtain are new and exhibit explicit invariance under 



renormalization group transformations. In particular, we include a correction of order 
to the fermion mass which has been left out in some previous analysis, but which is needed to 
ensure the explicit renormalization group invariance of the effective potential. In Sect. 
we extend the calculation of the effective potential to finite temperature. We analyze in 
particular the cancellation of ultraviolet divergences, and verify that indeed these cancella- 
tions take place without the need for counterterms other than the zero temperature ones. 
In Sect. we present a physical discussion of the thermodynamical properties of the sys- 
tem. We first analyze the mean field, or large A^, approximation and recover known results 
concerning chiral symmetry breaking, and its restoration at a finite temperature Tc. Then 
we show how the expansion breaks down for temperatures below T^. This is signaled 
in particular by the appearance of a Landau pole which invalidates the calculation of the 
effective potential for small values of the field. Finally we turn to the high temperature 
regime, where we find that the expansion becomes again a regular expansion. There 
chiral symmetry is restored and no Landau pole occurs. In this regime, the system behaves 
as a system of weakly interacting massless fermions, the corrections providing interac- 
tions which decrease logarithmically with increasing temperature. This is as expected from 
asymptotic freedom. However, since the coupling constant does not appear as an explicit 
parameter in our expression of the effective potential, this result is obtained only after a 
detailed analysis of the high temperature behavior of the effective potential obtained in the 
1/N expansion, and a comparison with the first orders of ordinary perturbation theory. We 
show that both approaches yield identical results when the coupling is the running coupling 
at a scale of the order of the temperature. 

II. THE GROSS-NEVEU MODEL. GENERALITIES 

The lagrangian of the Gross- Neveu model Jll 

£(^,V^)=^^^^ + ^(V^V^)^ (1) 
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describes interacting massless fermions in one spatial dimension. The summation over 
the flavors is imphcit in Eq. (|I]), e.g. ^if) = J2a=i4'a'4'a- As usual, ^ = 'j^d^, where the 
7 matrices are 2x2 matrices satisfying ■j^'ju + ■Julfi = '^dfiu, with g^j^^, = diag(l,— 1), and 
75 = 75 = 7°7^- 

The Lagrangian (|I]) is invariant under the discrete chiral transformation 

^p 757/^, i> -V'75- (2) 

Since tpip —ipip under this transformation, while the Lagrangian (|l|) remains invariant, 
the quark condensate (ipip) plays the role of an order parameter: its non vanishing indicates 
spontaneous chiral symmetry breaking, a situation met in the vacuum state 111]. 



In order to study the thermodynamics, we use the imaginary time formalism ^ 
and write the partition function Z as the following path integral: 

r/3 



Z = A/""^ J VipVip exp I - ^ dr j dx 



(3) 



where /3 = 1/T is the inverse of the temperature. We shall sometimes denote the space-time 
volume by / d^x = with L the length of the system (eventually to be taken infinite). 
In Eq. @ /iq = — 7°7^c?x is the hamiltonian of a free Dirac particle. The fields ip{T,x) and 
■iP{t,x) are antiperiodic, with period (3: tp{l3,x) = —tp{0,x). The (infinite) normalization 
constant M in Eq. (H), which depends on the temperature but not on the coupling constant, 
can be eliminated when necessary (see e.g. Eq. (^ below) by dividing Z by the (known) 
partition function Zq of free fermions at the same temperature: 



with Cp = \p\, and the factor 2 in the first term accounts for quarks and antiquarks. The 
last term (in parenthesis) is infinite, but independent of the temperature; it contributes only 
to the vacuum energy, and can be discarded. This is a trivial part of the renormalization 
which will be discussed at length in Sect. |T|. 

At this point, let us briefly digress on the notation that we shall use throughout to evaluate 
momentum integrals in the imaginary time formalism. Such integrals involve actually a sum 
over Matsubara frequencies and a true integral over the one-dimensional momentum p. They 
will be denoted by: 

n,odd-^^°° '^^ '' n, e?;en 27r 



where P denotes a 2-momentum, and the notation Y.n,odd (Jin, even) indicates a sum over 
fermionic (bosonic) Matsubara frequencies: u;„ = {2n + l)7rT for fermions, u;„ = 2mrT 
for bosons. Depending upon the context, the 2-momentum P will be either a Minkowski 
momentum or an Euclidean one. Whenever ambiguities may arise we shall use the more 
explicit notation Pm = for a Minkowski momentum, and Pe = {po,p) for an Euclidean 

one; we have = u"^ — p"^, and P| = Po + p"^- Functions of the 2-momentum P will be 
considered in general as functions of the Minkowski variables uj,p and, with a slight abuse in 
notation, denoted by either f{P) or f{u!,p); when Euclidean variables are appropriate, we 
shall write /{Pe) or f{ipo,p)- In doing finite temperature calculations one is led to set u = 
iuOni where uOn is a Matsubara frequency. Zero temperature contributions may be obtained 
from a finite temperature calculation by replacing TX^n by p^^doj / {2Tii) = /^dpo/(2vr), 
leading to an Euclidean integral / d^p/(27r)^. The way to handle the ultraviolet divergences 



will be described in Sect. pT[ in the present section we assume that all integrals are properly 
regularized, without making the regularization explicit. 



A. Effective potential and the expansion 

We now return to the partition function Z of Eq. (^). Following a standard procedure 
HI], we introduce an auxiliary scalar field a and write: 

Z=Af-'Jvae-'^I^''^'''z„, (6) 

where A/'ct = /r'cre"^/'^^^ <^^^ and the integration runs over periodic fields: a{P,x) = a{0,x). 
The quantity Z^^: 

Z^ = Af-^ J D7/>P^/>e-/'^'"^'(^-+'^('^))'^, (7) 

where h{a) = ho + -^1o, may be viewed as the partition function for a system of massless 
fermions in an "external field" a. We denote by 6*0- the fermion propagator in this external 
field; it obeys the equation 

{dri + h{ai)) S^{ti,Xi;t2,X2) = 5{ti - T2)5{xi - X2), (8) 

where ai = (j{ti,xi). For a given field a, the Gaussian integral in Eq. can be calculated 
in terms of 5*^. By taking the ratio Z„j Zq one eliminates the infinite normalization constant 
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and gets: 

In(^) =Trln5;i-Trln5o-\ (9) 

where the propagator 5*0 satisfies Eq. (H) with g = 0, and the symbol Tr imphes a trace over 
the Dirac matrices and an integration over space-time coordinates. 

The discrete chiral symmetry of the massless Dirac hamiltonian, Eq. (0), entails the 
following property: 75/i(cr)75 = h{—a), from which it follows that Z^^ is invariant under the 
transformation cr — > — cr. Since the weight function in Eq. (^) is also an even function of a, 
one expects the average value of a to vanish, unless there is spontaneous symmetry breaking. 

Whether symmetry breaking occurs or not can be deduced from the effective potential 
for the field a. To get this effective potential, we evaluate first the partition function in the 
presence of an external source j coupled to a: 

Z[j] = e-'^t^l = Af-' I PaZ,e-/^'"(5-'+i-). (10) 

The expectation value of a in equilibrium, that is, in the state of the system corresponding 
to the minimum of the free energy in the presence of the source j, is given by: 

A Legendre transform allows us to eliminate the source j in favor of the expectation value 
of the field: 

T[a^] = W[j]- j d^xa,{x)3{x). (12) 

When (jj is constant in space and time, we define (dropping now the subscript j): 

T[a] = V{a) J d'^x = (3LV{a), (13) 

where V^(cx) is the effective potential. Note that while F is dimensionless, V has the dimension 
of an energy density: it is the free energy per unit length (or minus the pressure) for a 
prescribed value of a. Note also that the discrete chiral symmetry implies that V is an even 
function of a, i.e., V{a) = V{—cr). 

The effective potential allows a simple determination of the equilibrium state. Indeed, 
since by construction 6T/6a{x) = — j(x), the equilibrium state in the absence of the source 
(i.e, for j = 0) is determined by the equation 



Furthermore, since the effective potential is a convex function (this follows from the general 
properties of the Legendre transform and the convexity of the solutions to Eq. (p!^ 

correspond to minima of the effective potential. In the absence of symmetry breaking, 
V{a) has a unique minimum, a = 0. When spontaneous symmetry breaking occurs, two 
degenerate minima appear, at nonzero values of a, say ±amin- 

Strictly speaking, the technique of Legendre transforms allows us to construct V"((t) as 
indicated above for all values of a outside the interval [—amin, +o'min], and yields a constant 
value V{a) = V{amin) inside this interval (for a pedagogical discussion of this point see e.g. 
|3^). Note that this construction is somewhat formal; in particular the values of a in the 
interval [—cTmin, +(^min] are not reached as minimal free energy solutions in the presence of 
a constant source j. Now, in the approximations to be developed shortly, we shall find that 
the relation between the source j and a is multivalued: one solution corresponds to the 
minimal free energy for a given j, while other solutions yield higher free energies and values 
of a inside the interval [—(Tmin, +o'min]- By keeping all solutions, we define a continuation 
of the effective potential which is not a constant in the interval [—amin, +o'min]- It is this 
continuation of the effective potential that is used in particular to calculate the fluctuations 
of the field a. 

In this paper we consider the first two orders in the expansion of the effective 
potential V. These are obtained by evaluating the path integral (|TDp in a saddle point 
approximation |^|, and then performing a Legendre transform. The leading contribution 
will come from the saddle point itself, the correction from integrating the fluctuations 
about the saddle point. 

The value of the fleld a at the saddle point, i.e., the value cXc for which the exponent in 
Eq. (0) (with Zcr of Eq. (|^)) is extremum, is the solution of the following self-consistent 
equation (commonly referred to as the "gap equation"): 

1 dlnZ, 



PL da 



^ m.., (15) 



where we have used Eq. (0) to express the derivative din Zcr/da\^^ in terms of the quark 
condensate {ipip)^^^, and we restrict ourselves to solutions which are constant in space and 
time. The quark condensate {ip'ip)a^ is easily calculated with the help of the propagator 
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(see Eq. (P)), whose Fourier transform for constant a reads: 



This is the propagator of a quark with mass M defined as: 



M = 4^. (17) 



For a given (constant) value of cr, the quark condensates reads then: 

1 



(V'V'). = - / {d'P}Ti 



-f + M 

from which we get {ipil))^^ by setting a = (Tc- From the equations above, we observe that 
(V'V')crc is of order A^, so that is of order \/N (the source j is to be considered as a constant 
of order y/N); as for the mass parameter M, it is independent of A^. 

At this level of approximation, the free energy in the presence of the source, W[j] = 
— ln^[j], is obtained from the value at the saddle point of the integrand in Eq. (p^O]): 



W^%] = {^\a!+ja^)f3L-\n^, (19) 

and, using Eq. (0), the expectation of the field is simply a = ac- Note that the gap equation 
(P^) may have multiple constant solutions corresponding to the same constant j (this can 
be seen explicitly from the expressions given in subsection |III A| ). Keeping, as discussed 
above, all solutions (i.e., not only that with the lowest free energy), and eliminating the 
source j according to Eq. ([121) , one obtains the following "continuation" of the effective 
potential: 

As we shall see in subsection pi A| , this function has two degenerate minima, but is not flat 
between the two minima: the values of a within the two minima correspond to the solutions 
of the gap equation which, for a given value of j, accompany the solutions with minimal free 
energy. 

In order to get the next order in the 1/A^ expansion, we expand the field a around the 
value (Tc, i.e., we set a = o'c + d', and do the corresponding change of variable in the functional 
integral (0). The terms linear in a cancel because ac satisfies Eq. ([15|). There remains then 
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FIG. 1: The contributions to the free energy in the expansion: the left diagram represents 
the Hartree contribution; the right diagram is one of the infinite set of "ring diagrams" . 

in (p!OD a gaussian integral over the fluctuations of the field a. This is easily evaluated, with 
the result: 

-^^ j Va exp J [d'Q] a{Q) (l + g'U{Q)) a(-g)} 

= exp|-^ l[d'Q]\n{l+g'niQ))y (21) 

The quantity Il{Q) that appears in this equation is the contribution of the one loop diagram: 

n(Q) = l/KA-}Tv[^^^^]. (22) 

which is calculated in App. A; it plays the role of a self-energy for the a meson, whose 
propagator D{Q) can be read from Eq. (|2TD above: 



Combining the value W^^^ obtained above, Eq. ([19|) , with the gaussian integral (|2l|) , we 
obtain W[j] at order 1/A^: 

W[J] = {l^l + J^c) - 1^ ^ + ^ / id'Q] In (l + g'm)) = W^'^ + W^'l (24) 

By taking the derivative of W[j] in Eq. (|2^) with respect to j, one obtains a, the expectation 
value of the field a. (Note that cTc depends implicitly on j through the gap equation ([T5|).) We 
can write a = (Tc + Sa, and it follows immediately from Eqs. (p^ and ([T5| ) that 6a is of order 
relative to cXc. In computing the Legendre transform T[a] = W^'^^ + W^^^ — ja(3L (see 
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Eq. (|T2D), we may simply replace cXc by a, since the error made is of order (the terms 

linear in 5a drop in because Uc satisfies the gap equation (0); similarly, one makes an 
error of order l/N"^ in ignoring 5a in W^^"* which is already a quantity of order l/N). The 
elimination of j is then straightforward and we finally obtain the effective potential in the 
form: 

V{M) = N^-nJ {d'P} In (l - ^) + ^ / [d'Q] In (l + /n(g)) , (25) 
where M is given by Eq. ( P^Tj ) with a replaced by a . We shall write this effective potential 



as 



V{M) = NV^^\M) + 1/(1) (M), (26) 

where NV^^\M) is the leading order contribution and V^^\M) the next-to-leading order 
one. Before renormalization, NV^^\M) is the sum of the first two terms in Eq. (|25|) , and 
V'^'^\M) is the last term. After renormalization, as we shall see, the first term in Eq. (|25D 
contributes also at order For this reason, we shall often refer to the first two terms 

in Eq. (^) as the "fermionic" contribution, and to the last term as the "bosonic" one: the 
first two terms represent indeed the free energy density of massive fermions in the "Hartree 
approximation" |3y] , while the last term may be viewed as a the one-loop contribution of the 



bosonic degrees of freedom associated with the a field. In terms of the fermionic variables 
of the original lagrangian, the corresponding contributions to the free energy can be given a 
simple diagrammatic interpretation: the first diagram in Fig. corresponds to the Hartree 
contribution, the second is one of the family of "ring diagrams" representing the bosonic 
contribution. 

Before closing this subsection, let us recall that the effective potential is the generating 
functional of the irreducible n-point functions for the a field at zero momentum. In partic- 
ular, the following, leading order, relation for the 2-point function can be easily established 
at zero temperature: 



d^l/W 



D-\Q = Q-M) (27) 



dM2 

where = 1 + g^Il is the inverse of the a propagator (see Eq. (|23|) ). The same relation 
holds at finite temperature, but care must be exerted in specifying the limit Q — > which 
is in general non analytic (it depends on which limit a; — > and g ^ is taken first). 
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(a) (b) (c) 

FIG. 2: The fermion condensate {ip^): (a) leading order; (b) correction, the dashed line 
represents the a propagator; (c) same as (b) in the fermionic language. 



Thus, at finite temperature, the right hand side of the relation (|27|) has to be understood as 
\imq^QD^^{uj = 0,q; M) (see subsection [VB| and App. A). 



B. Gap equation and quark condensate 

As we have argued before, the equihbrium state is obtained by minimizing the effective 
potential with respect to M. One gets then the gap equation in the form: 







dM 



= ^^ + (^^)m_, (28) 

M 9 



where we have used the fact that (ipip) may be obtained as the derivative with respect to 
M of the last two terms in Eq. ( [25| ) (this may be seen by going back to the expression (^ 
of the partition function and using the fact that (ipip) can be obtained by differentiating Za- 
in Eq. (|^ with respect to a; see also Eq. ([I5|) ). Thus, the value of M at the minimum of 
the potential and the value of the quark-condensate in equilibrium are proportional: 

M^,„ = -|^(V^^)m_, (29) 

which indicates that a (or M) and {ipip) cii'g equivalent order parameters for the discrete 
chiral symmetry. 

We shall solve Eq. ( pSj ) order by order in the 1/A^ expansion. To do so, we write: 

M™„ = Mil + ^Mii. (30) 



In leading order, we have: 







dM 



(0) N 
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(a) (b) (c) (d) 

FIG. 3: Contributions to the fermion mass, (a) leading order; (b) correction to the minimum 
of the effective potential; equivalently, the sum of the contributions (a) and (b) can be represented 
by diagram (c) where the insertion with a cross represents Mmin- Diagram (d) represents an 
additional correction coming from S. As in Fig. ^, the dashed line in (d) represents the a 
propagator 

where (■?/''?/') ^(o) is obtained from Eq. (|^) with M^j^ substituted for M. This gap equation 

min 

needs to be solved exactly in order to get M^j^^. The correction M^^}^ is obtained by solving 
Eq. ( p8D to accuracy. We get then: 



mm 







(IM 




d2y(o) 




dM2 





(32) 



To the expansion (|30D of Mmin corresponds, according to Eq. (l29|), a similar expansion of 
the quark condensate: 

{ijifj) = N{i)i))^^'^ + (33) 

The two contributions to the quark condensate correspond to the diagrams displayed in 
Fig. 0. The leading order {ipip)^^^ corresponds to the contribution of the first diagram 
evaluated for Mmin = ^min- The order correction [ifjil))^^^ contains the contribution 

- ' 'ye-Q] , '""'i'L (34) 



dM 2J ' 1+ g'^Yi{Q) 
associated with either of the last two diagrams in Fig. |^. In addition, as explicitly indicated 
in Eq. (|3^) , {ijj-ip)^^^ receives also a contribution from the first diagram in Fig. |] evaluated 



at Mmin ^ M, 



(0) 



mm' 



Before closing this subsection, it is worth emphasizing that neither the minimum Mmin 
of the effective potential, nor the quark condensate are physical observables: after renor- 
malization, and beyond leading order in the expansion, their values will depend on the 
renormalization scale. 
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C. The quark mass Mf 



As we have mentioned, chiral symmetry is spontaneously broken in the ground state, 
and because of their couphng to the condensate the quarks acquire a mass. This mass is 
a physical quantity that we shall denote throughout as Mj. That is, Mf will consistently 
refer to the mass of the quark in the vacuum, a quantity to be kept constant at all orders of 
our approximations. 

In leading order in the expansion, the coupling to the condensate is the only contri- 
bution to the fermion mass and we can identify M^}^ with Mf. This identification between 
Mmin and Mf does not hold in higher order. Starting at order there are other con- 

tributions to the fermion mass besides M^m- Quite generally, the fermion mass is given by 
the pole of the fermion propagator at P| = —Mj. The fermion propagator itself, S{P), can 
be written as 

S-\P) = S;\P) + S(P) = -r + M™„ + S(P). (35) 

The various Feynman diagrams contributing to the fermion mass are displayed in Fig. |^. In 
leading order, S(P) = 0. At next-to-leading order, S(P) can be written as: 

S(P) = af + b, (36) 

where a and b are functions of P given in App. C. The equation S^^{^ = Mf) = yields 
then Mf as a sum of two contributions: 

Mf = Mmin + Mj:, (37) 

where Ms (which is of order is calculated in App. C. 

At finite temperature, we shall define a temperature dependent mass Mf{T) by the 
equation: 

S-\uj = Mf{T),p = 0;T) =0. (38) 

where S{uj,p;T) is the fermion propagator at finite temperature. At leading order in the 
expansion Ti{uj,p;T) = 0, and one can simply identify the fermion mass Mf{T) with 
the minimum of the finite temperature effective potential. At next-to-leading order, the 
evaluation of the Mf{T) would require the calculation of T,{ijj,p] T). But, as we shall see in 
the following sections, this is not needed to evaluate the pressure at next-to-leading order. 
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III. THE RENORMALIZED EFFECTIVE POTENTIAL AT T = 



As we have already mentioned, the formulae given in the previous section contain ultra- 
violet divergent integrals. We discuss now the procedure of renormalization which allows us 
to get rid of the infinities. First we need to specify our regularization: this will consist in a 
simple cut-off A on the length of the Euclidean momenta (for a discussion of more sophisti- 
cated regularizations in this type of models, see e.g. Ref. We then obtain the effective 
potential as a function Vb{Mb', Qb, -^), referred to as the "bare" potential; the variables as 
(or Mb) and qb, are the bare field (or mass) and coupling constant, respectively. 

Divergent integrals appear in dV^/dM^ and d^Ve/dM^, i.e., in the one and two-point 
functions of the a-field. The corresponding diagrams at leading order and at next-to-leading 
order are displayed in Fig. |. Besides, there are divergences which are independent of Mb', 
these are eliminated by subtracting from Vb{Mb; gs, a constant C{gB,-^)- 

In order to construct the renormalized effective potential, we define a renormalized field 
a, a renormalized coupling constant g, and a renormalized mass M (related to a and g by 
Eq. (|13)). We set: 

(TB = ^(r, g^ = ^g^ Mb = ^'M. (39) 

where the renormalization constants Z and Z' are dimensionless functions of the renormal- 
ized coupling g, and the cut-off A. These functions may be expanded in powers of 

Z = + Z' = Z'W + Iz'(i), (40) 

and chosen so as to absorb the ultraviolet divergences at each order. Because they are 
dimensionless, Z and Z' can actually depend only on the ratio of A to another mass scale Mq, 
to be referred to as the renormalization scale, at which the renormalized n-point functions 
are specified (see below). In fact, we shall find convenient to consider Z and Z' as functions of 
both A/Mf, where Mf is the fermion mass, and M^/Mf (the latter quantity being eventually 
a finite function of the renormalized coupling): the dependence in A will be fixed by the 
elimination of ultraviolet divergences, that on Mq will be determined by renormalization 
conditions that we now present. 

To this aim, we note first that the correction to the vertex aipip is of order (see 
Fig. § f). In leading order, this vertex is therefore not renormalized and we may choose 
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(a) 



(b) 




(c) 





(d) 



(e) 




(f) 



FIG. 4: One-point (a,c) and two-point (b,d,e) functions for the a field, in leading order in the 
expansion (a,b), and at order (c,d,e). Diagram (f) is the first correction to the quark-sigma 
vertex; it is of order The dashed lines represent a propagators. 

to balance the renormalization of the field with that of the coupling constant, i.e., set 
gB<^B = or equivalently, according to Eq. (|39D , set Z'^^'^ = 1. As for the constant Z^^\ it 
will be determined by the following condition: 

d2\/(0) 



dM2 



Mo ^ 

According to the relation ( pTf ) (extended to the corresponding renormalized quantities) this is 
equivalent to the condition D{Q^ = 0; Mq) = 1 for the renormalized a propagator. At next- 
to-leading order we shall impose |^ that the first and second derivatives of the potential 
at M = Mo are not changed by the corrections, that is: 



dV(i) 



dM2 



0, 



d2V(l) 



Mo 



dM2 



0. 



(42) 



Mo 



Summarizing the procedure that we have outlined, we may write the renormalized effec- 
tive potential V as follows: 



V{M; g, Mo) = Vb(Mb; Qb, A) - C{gB, A). 



(43) 



By construction, the right hand side has a finite limit when the cut-off A is sent to infinity. 
For large A, the renormalized potential becomes therefore independent of A, i.e., it satisfies 
the renormalization group equation AdV^/dA = 0, where the derivative is taken at fixed 
values of the renormalized quantities, and fixed Mo. This can be written explicitly as 

d 



i{9b)M^ 



B 



dM, 



B 



{VB{MB;9B,A)-C{gB,A)) = 0, 



(44) 
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where 

dgB QB d\n{Z'/Z) A OMb I dlnZ' 

Alternatively, expressing the fact that the renormalized potential in Eq. ( ^31) does not 
depend on the scale Mq, one can write the renormalization group equation as Mq dV"/dMo = 
0, where the derivative is now taken at fixed bare quantities, and fixed cut-off A. Explicitly: 



ViM;g,Mo)=0, (46) 



where 



dg g d\n{Z/Z') M, DM 1 dlr^Z' 

These various functions will be calculated in subsection |III C| , where we shall also dis- 
cuss some consequences of the renormalization group equations. We note here a general 
feature of the model that is useful to keep in mind in order to understand the logic of the 
construction of the effective potential in the next subsections. The model depends initially 
on two parameters, the bare coupling gB and the cut-off A. Since these will be adjusted 
so as to reproduce the fermion mass M/, gs will become effectively a function of A. This 
relation between gB and A depends on the accuracy with which Mf is calculated, and is 
therefore modified at each order in the expansion. The same property holds in the 
renormalized theory: there exists a relation between the renormalized coupling g and the 
scale Mo, which is redefined at each order so that the calculation of the fermion mass at 
that order reproduces the value Mj. 

We shall, in the next two subsections, construct the leading and next-to-leading order 
contributions to the renormalized effective potential. We shall see that it is possible to write 
these in terms of a variable Xm related to M, but which is invariant under renormalization 
group transformations (M is not invariant beyond leading order). This choice of variable 
will make it obvious that the potential is independent of Mq, and it follows from the relation 
between Mq and g alluded to above, that it will be also independent of the renormalized 
coupling. 
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A. The leading order contribution 



The zeroth order in the expansion of the potential, is the fermionic contribution 
in Eq. (pS]). At T=0 it can be written as: 



NV^'\Mb) = N^Ml-N J 



In 1 - 



p2 



(48) 



This expression is divergent. Introducing a cut-off A on the length of the Euclidean momen- 
tum < A), one obtains: 



<'(A/i,) 



1b 



O ' ^ 



2gl ■ Atx A2 ^/ ' " V A2 ' ■ ^'^^^ 
We proceed now to the renormalization of the potential along the lines described in the 
previous section. Since at leading order Z' = Z'^^^ = l (see the discussion after Eq. (^OD), we 
can simply replace Mb by M, and we need only to renormalize the coupling constant. After 
this, there will be no need for further subtraction. By expressing (^s in terms of g in Eq. 
(H), using Eq. (|3|) with Z = Z^^\ we get: 



M = — + In — - 1 

^ ' 2c/2 47r I A2 



(50) 



The A-dependent term is eliminated by choosing 

, A2 



where the finite constant Z'^'^'^ is fixed by the renormalization condition (p|): 

Z(o)(Mo,^7) = 1 -//7r + (//7r)ln(M;/Mo). (52) 



The renormalized potential at leading order can then be written as: 



M2 



M2 


/ M2 \ 


M2 / 




n ^ - 1 




47r 
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It has two minima, at ±M^}^, with: 



M;^l{Mo,g) 



M2 
Mn 



Mne^-'/s'. 



3 + ^ 



27^ 



9 



(53) 



(54) 



The existence of such non trivial minima of the effective potential indicates spontaneous 
symmetry breaking with, according to Eq. (p9D , a non vanishing value of the quark conden- 
sate. 
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It remains to verify that V^^^ does not depend on Mq. To do so, we remind first that, at 
this order, we can set M^^}^ = Mf, or equivalently (see Eq. (||)) Z^o) = 0. Then, one gets: 

V^'\M) = J^[^^M]-^)= j^XMilnxM - 1), (55) 

where we have set 

xm = My Ml (56) 

Since at this order M is not renormahzed, the variable xm is clearly a renormalization group 
invariant, and so is the effective potential. At next-to-leading order, we shall see that it is 
still possible to express the effective potential in terms of a renormalization group invariant 
Xm whose definition generalizes that in Eq. (^) so as to include an Mo-dependent correction 
of order (this correction will compensate the Mo-dependence of M; see Eq. (|5B|)). We 
note also that, as expected from the discussion at the end of the previous subsection, the 



expression (|55|) of the effective potential does not depend on the coupling constant anymore: 
by eliminating the explicit dependence on Mo, we have also eliminated the dependence on 
9- 

In order to get the expression ( ^5]) of the effective potential, we did not need the explicit 
expression of given by Eq. (|52D. This is needed however to specify the explicit relation 
between Mo and g. From the equation Mj = M^}^{Mo, g), or equivalently Z^^^{Mo,g) = 0, 
one gets: 

«^<^« = rTMlVM^- 

This formula indicates that g becomes vanishingly small as Mo ^ oo, as expected from 
asymptotic freedom. It also shows that g diverges when Mo Mj/e, a property that we 
shall discuss further shortly. 

Before we do that, we note that the main results obtained so far in this subsection could 
be derived by working solely with the bare quantities. In terms of the bare parameters and 
the cut-off A, the minima of vI^\Mb) are given by the solutions of the gap equation: 

Mr Mr / A^ ^ 



B 

gl A^"'\Mj, 



2 - • (58) 



The right hand side of this equation is —{il)ip)/N, i.e., it is proportional to the quark con- 
densate (see Eq. (p8D). Aside from the trivial solution Mb = 0, the gap equation has two 
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FIG. 5: The quark-quark scattering amplitude corresponding to the exchange of a a meson (left- 
hand side), and its expression in terms of the fermionic variables. 

degenerate solutions at Mb = ±M^j-„ with: 

Ml^l{A,gB)=Ae-^/^B. (59) 

By identifying one of the non trivial solutions with the fermion mass, as we did above, e.g. 
setting Mf = Mj^^^{A, gs), we get the relation between the bare coupling and the cut-off: 

This may be used to eliminate A from the expression (|49| ) of the potential and recover 
Eq. (H) (recall that at this order Mb = M). 

We end this subsection by considering the quark-quark scattering amplitude T{uj, g; M), 
where the quark mass M is considered as an independent parameter. In leading order (in 
perturbation theory) this scattering amplitude is simply the bare coupling constant squared 
Qb- Including the exchange of the a meson, which contributes at the same order as the bare 
coupling in the expansion, one obtains: 

A diagrammatic interpretation of T{Qe',M) is given in Fig. |^. The scattering amplitude 
can also be expressed in terms of renormalized quantities (see App. A): 

r{QE;M)=g'Do{QE;M), (62) 

where is the square of the renormalized coupling constant at the scale Mq, i.e., g^ = 
^'^(Mq) (given by Eq. (^7|)) and Dq is the renormalized a propagator (given by Eq. (|A6|) ). 
Note that the explicit g"^ in the r.h.s. of Eq. (|6^) cancels with the factor l/g"^ contained in 
Dq (see Eq. (|X6|)) so that T{Qe', M) is independent of the renormalisation scale Mq, i.e., the 
scattering amplitude is a renormalization group invariant. Using the explicit expression of 
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Do{Qe', M) given in App. A, one easily shows that T{Qe = 0; Mf) = g'^{M^), and that at 
large Ql (Q| > M|, M < Mf), r{QE\ Mf) ^ 27r/ln(Q|/M|) ^ g\Ql). This logarithmic 
decrease of the interaction strength at large Q\ reflects of course the property of asymptotic 
freedom. 

As explained in App. A (see Eq. ( |A^ ) and the discussion that follows), the a propagator 
develops a pole at flnite Qe when M < Mf /e and, according to Eq. (§^), so does the quark- 
quark scattering amplitude. This pole which does not correspond to a physical excitation 
of the system is commonly referred to as a Landau pole; as we shall see, it is responsible 
for several difficulties in our next-to-leading order calculation at flnite temperature. We 
shall denote by M^, the value of M at which the Landau pole appears at = 0. Here 
M^ = Mf/e. As we have seen before, at = the scattering amplitude is simply T{Qe = 
0;M) = ^2(M2) which diverges when M ^ M^ (see Eq. (P)). Since T{Qe = 0; M) 
is nothing but the inverse of the second derivative of the leading order effective potential 
V"(°) (see Eqs. {^) and (H), and also (^)), M, is also the zero of the second derivative of 
(^). At zero temperature, M^, is far from the physical quark mass Mf, and the quantum 
fluctuations of the a fleld do not probe values of M close to M^,; the Landau pole is then 
harmless. We shall see that this is no longer the case at flnite temperature. 



B. The 1/N contribution 



are two such contributions that we call vj^^ and Vf, 



We consider now the next-to-leading order contributions to the effective potential. There 

The latter is the bosonic contribution, 

i.e., the last term in Eq. (pS)). In terms of bare quantities, this reads: 

1 cPq 
2J (2^ 



V^^,{MB■.gB,^) 



In 



l + f7|no(QE;MB,A) 



(63) 



where the cut-off A applies to the magnitude \jQ% of the Euclidean momentum, and 
^f){QE'i Mb, A) is given in App. A. The other contribution, Vf^\ originates from the renor- 
malization of the coupling constant qb and the mass Mb, at next-to-leading order, in the 
fermionic contribution, Eq. (^9]) . It will be examined later in this section. 

We flrst address the problem of flnding a constant C{gB,-A) to eliminate the M- 
independent divergences ofV^^)^. We deflne: 



C{gB,A) 



In 



l + glUo{QE;M^,A) 



Ml 
An 



In 



A^ 
Ml 



Ml 
An 



In In ■ 



A 

Ma' 



(64) 
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with Ma = Ae~'^/^s^ and we verify now that this subtraction fulfills our requirements. First, 
we express vI;^^{Mb'-, Qb-, ^) — C{gB,A) in terms of the renormalized mass and coupling 
constant according to Eqs. (|39D. Note that since V^^^ is already of order we need only 
the leading order renormalization constants. That is, we can set Mb = M and = g'^/Z^^\ 
with given by Eq. ( pi]) with Z^^^ = 0. Furthermore, we can also replace Ma = Ae '^^^b 
by Mf, using the relation (^) and the identification at leading order of M^]„ with Mf. 
Thus, at this stage: 



'> ^ ' 2 J 27r 2 



D^\Qe;M) 



Mj Mj A 
H ^In — 2 ^Inln . (65) 



Do\QE;Mf)\ ' 47r"^M2 ' Att """ Mf 



{2ny 

where Dq\Qe; M) = + g^Uo{QE; M, A) is the renormalized inverse a propagator given 
by Eq. (^. 

The integrand in Eq. ( |65D vanishes when Q^^ oo, but not fast enough to make the 
integral convergent as A — > oo. Indeed, from (|A8|) one gets 



where: 



Do\Qe;M) ^ 2A{M) 2B{M) / M^ \ 

D,\QE;Mf)- ^ Ql ^ Ql\n{Ql/Mf)^''[Q%ln{Ql/Mj))' ^^^^ 



A{M) = M' -Mf, B{M) = M' - Mf - M' \n{MyMf) . (67) 

The terms proportional to A and B yield the following contributions to the integral: 

A(M)^ A2 S(M), , A 

^^In— ^, ^-^Inln— , (68) 

47r M| Att Mf ^ ^ 

which contain M-dependent but also M-independent divergences. The latter cancel exactly 
the last two terms of Eq. (|65D. The former contribute M-dependent terms to Eq. (|65D 
which are of the same form as those of the leading order potential, V^^\ namely M^ and 
M^lnM^ (see Eq. (0)), and they can be eliminated by an appropriate adjustment of the 
two renormalization constants Z and Z'. To see that, we need to return to the leading order 
contribution to the bare potential, (Eq. (^)), and trade then Mb and gB for M and g, 
using Eqs. (p9|). We get: 

By expanding Z and Z' up to order 1 /A^ as in Eq. (^Oj) , one can write the above expression 
(H) as + {\lN)vf\ with given by Eq. (|5|) and, dropping terms of order 
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As mentioned at the beginning of this section, V^^\m) (from the equation above) adds 
to V^^\m) (from Eq. (|65|)) to yield the complete next-to- leading order contribution to the 
renormalized potential, V'^'^^M) = vj^\M) + V^^\m). The requirement that A disappears 
from V"*-^-* (M) then allows us to determine the A-dependence of the renormalization constants 
and Z'^^l One gets: 

,2 



2tx 



(71) 



Z'«=lnln(^A^+Z'«, (72) 

where Z'^^^ and Z^^^ are finite constants to be fixed using the renormalization conditions 
(^). We postpone this determination till later. Keeping along Z''^^^ and one can then 
write the renormalized effective potential in the form: 

,.a.(M) . z<..|! + z...^^ ,„ |! _ _L(M^ - MJ) + f f„(.„). (73) 

with xm = M'^/Mj (see Eq. (^)), and Fq{x) is a finite function defined in App. B. 

The renormalization conditions (^21) needed to fix the finite constants Z^^^ and Z'^^^ 



require two derivatives of V^*^^-*. The first derivative of V"*-^-* with respect to is easily 
obtained from ([73|) : 

(^V^(i) £(1) 1 1 
where Fq = dF^/dx. The second derivative of V"*^^) may be calculated as 

One thus easily gets: 

^2t/(1) 7/(1) 11 

-TXTT = ^ + ^(In^M + 3) - - + -F^(xm) + ^K(xm), (76) 

where F^' = d^Fo/dx^. The derivatives of Fq{x) are given in App. B. 

In this subsection we shall in fact only fix Z^^^ and leave Z'^^^ undetermined. The reason 
is that, as we shall see, it is not necessary to specify the value of Z'^^^ in order to get 
a renormalization group invariant expression for the effective potential. By imposing the 
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second renormalization condition (|42|) , viz d?V'^^^ / dM'^ = at M = Mq, we can express 
in terms of Z''^^^ Using Eq. ([76|) , one gets: 



27r 



2 _ (3 + lna;o) - F^(xo) - 2xoK(xo) 



where the quantity Xq is the following function of g: 



(77) 



(78) 



At leading order xq = M^/Mj, as can be easily deduced from Eq. (^7|). We can then write 
the next-to-leading order contribution to the renormalized potential as follows: 



y«(M) 



M 



Z'(iM In — 2 - 3 - In xo - {xo ) - 2a;oF^' (xq) 



M2 M2 



The full renormalized potential up to next-to-leading order is thus V = NV^^^ + V^^^ with 
V"(°) and 1/(1) given by Eqs. (^) and (|79D , respectively. 

The effective potential that we just constructed has an apparent dependence on the scale 
Mq. However, as we did in the leading order case, it is possible to exploit the relation 
between the fermion mass Mf and the value Mmin of M at the minimum of the effective 
potential, in order to get rid of this apparent Mq dependence. The situation is somewhat 
more complicated here because, as explained in subsection |IIB| , the fermion mass receives 
two contributions at order Mmin and M^. The latter is calculated in the App. C and 
takes the form: 



Ms = M 



f 



2N N . 



(80) 



where cp is a numerical constant, = —0.126229 (see Eq. ( pi6| )) . The value of Mmin can 
be calculated by setting M™„ = M^^^l + {l/N)Mi^l (see Eq. (|0|)), with M^^l and M^^l 
given by Eqs. (Q) and (|32D, respectively. The required derivative of V^^^ can be obtained 
from Eq. (^ (together with (jT^)). Note also that M^]„ = Mj + 0{1/N) and that, at the 
mmimum, d2r(o)/dM2 = l/vr. One thus gets: 



M. 



(1) 



(0) 



^ Inxo + l) Z'^'^ + iF^(xo) + xoF^xo) - ^F^(l) 



(81) 



Remembering that Mf = Mmin + Mj^, (see Eq. (0)), and combining Eqs. (pO| ) and (0), one 
can therefore write Mf as: 



Mf = Mo 



L 



1 + 



1 

N 



^(Inxo + 3) + ^F^(xo) + XoFaxo)-^F^(l)+y^ 



(82) 
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For the purpose of eliminating Mq, it is convenient to rewrite Eq. ( p2D as follows (neglecting 
terms of order l/N"^): 



^ Mf 1, 1 

In — - = — In xn H 

Mo 2 " I 2 



(lna;o + 3) + ^F^(xo) + xoF^Kxo) - ^F^(l) + • (83) 



We shall return to this equation in subsection |III C| below. For the moment we note that it 
allows us to easily eliminate the term ~ In Mq in the leading order potential V^*-'^-' (M) given 
in Eq. (|53|) . Combining the resulting expression with Eq. ([79D and droping again terms of 
order (in particular, using xq = M^/Mj + 0{1/N)), one finally gets 



V{M) 



NV^o) ^ ^(1) 

-1 + ln 



47r 



M2 
Ml 



+ 



2tx An 



Fo(xAf), (84) 



where ^ = 2^- F^(l) 



0.8661636. At this point we observe that 



1 



0{1/N^). (85) 



This allows us to express the full effective potential V entirely in terms of the redefined 
variable xm- 



Xm 



M2 



1 + 



aT 



(86) 



The final result for the effective potential can then be put in the form (with a slight abuse 
in the notation): 



M2 

V{xm) = -r- [Nxui^'^XM - 1) + i{xM - 1) + Fo{xm)] ■ 

4:71 



(87) 



where we have explicitly subtracted the constant 2 + ^ so that V^^\xm = 1) = 0. Note that 
the definition of xm in Eq- (§5) differs from that in Eq. (^) by a term of order 1/A^. There 
is no contradiction however since Fq{x) is already of order 1/A^, so that using xm = M'^/Mj 
instead of the definition (^) inside Fo{xm) is equivalent to ignoring terms of order l/N"^. We 
shall see in the next subsection that xm is a renormalization group invariant. The effective 
potential V{xm) is thus explicitly independent of the scale Mq, and for the reasons discussed 
at the beginning of this section (after Eq. (|47|)), also independent of g. Note that in terms of 
the redefined variable xm (Eq. (pB]) ), the leading order contribution in Eq. ( ^71) is identical 
to that of Eq. (|5D. 
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2 



N=100 




FIG. 6: The effective potential divided by NMj/Air as a function of xm for the values = 1, 
= 3, = 10 and N = 100. The minimum occurs at Xmin = 2.721 for = 1, Xmin = 1.145 for 
AT = 3, Xmin = 1-068 for A^ = 5, Xmin = 1-029 for A^ = 10, and Xmin ~ 1 - 2ip/N « 1 + 0.25/Ar 
for larger values of A^. The 1/A^ correction has been adjusted so as to vanish at xm = 1; it also 
changes sign at this point. 

The effective potential (^) is displayed in Fig. ^ for various values of A^. One sees that 
the correction becomes rapidly a small correction. For = 1, the approximation is 
of course meaningless, but already for iV = 3, the shape of the potential is qualitatively 
the same as at leading order. The variation of the minimum with may be obtained as 
explained in Sect, y from the general equation (p2D: 



X 



(1) 

min 



1 I fll'l'll- \ 







dxM 


2;m=i 


d2y{o) 






XM = l 



(88) 



From the exphcit expression ( ^7|) we get d'^V^'^y dxj^ ^ = Mj/An, and dV^^^dzu ^ = 
(M|/47r)(^ + F^(l)) = (pMj /2tt, so that, at next-to-leading order, x^in = 1 - '^f/N ^ 
1 -|- 0.25/A^. The difference between this value and that of the "exact" minimum found 
numerically from Eq. (|87D can be attributed to terms of order l/N"^. As can be deduced 
from the caption of Fig. ^ for A^ ^ 10 these are small, confirming the previous study of 
Ref. [Q. (Of course, such considerations say very little about the magnitude of the true 
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corrections.) It is interesting to recover the previous result for Xmin in a different way. 
Using the definition of xm in Eq- (|86|), we can write Mmin = \/Xmin Mf{l — Z'^^"* /2N). The 
equation Mf = Mmin + M^, reads then: 



Mf = M 



f 



1 + 



X, 



(1) 



2N 2N 



+ Mf 



2N ^ N 



(89) 



where we have used Eq. (^Up for Ms and the first equation in (^Sp. Eq. (^) gives back our 
previous resuh x^m = — 2(y9. It also clearly illustrates how the scale dependence cancels out 
between Mmin and Ms, whereas the magnitude of each individual contribution depends on 
Mo (through Z'^^)). We come back on this issue in the next subsection. 

Note finally that the potential is plotted only for xm > = xm,- This is because for 



Xm < a; A/* a Landau pole appears in the a propagator (see subsection |III A| ), making the 
function Fq complex (see Eq. (^) or, equivalently Eq. (pT])). The influence of this Landau 
pole on the shape of the effective potential is clearly visible in Fig. |]on the plot of V{xm) 
for = 1: not only does the potential becomes complex for xm < xm,, but the "1/A^ 
correction" is large and negative for xm ~ xm^- As soon as > 3, however, this peculiar 
behavior is no longer visible. 



C. Renormalisation group analysis 

In this subsection, we shall clarify the behavior under renormalization group transfor- 
mations of various quantities that have been introduced in the previous subsections. In 
particular we shall verify that the expression (|82|) of the fermion mass is invariant, and so is 
the variable xm defined in Eq. (0). But before we do that, let us verify that our calculation 
reproduces known results for the renormalization group functions. 

The functions /3 and 7 are most easily obtained by differentiating the renormalization 
constants with respect to the cut-off A, according to Eq. (^). Using Eqs. (0), ( [7T| ) and 
(^) one easily gets: 

P{9B) = P^'\gB) + j^P^'\gB), (90) 



with 
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Eqs. ( |90D and (0) show that, for = 1, the leading term in a small (^-expansion (i.e., 
that of order g^) in the /5-function vanishes. This is in agreement with earlier calculations 



and the fact that for = 1 the Gross-Neveu model reduces to the Thirring model |^ 
for which the /3-function vanishes identically. (The term of order does not vanish when 
= 1 in Eq. (^Tj) because the correction does not contain all the terms of order g^ in 



the small g limit). The two- loop /^-function has been calculated by Wetzel [B^ and reads 
f3{g) = -(1 - l/iV)(^727r) + (1/A^ - 1/A^2) ^^5/4^2))^ ^^^^ g^^^Uy ^^lat P^^^gs) follows 
immediately from the explicit relation (|60|) between gs and the cut-off A. 

Following similar steps, one obtains the function 7 = 'y^^^ + (1/A^)7^^''. Since at leading 
order, Z' = 1, 'j^^^gs) = 0. Using Eq. ( [72D and also the relation (^) between the cut-off A 
and the bare coupling constant we get: 



9l 



(92) 



It is instructive to repeat the calculation of the renormalization group functions by dif- 
ferentiating with respect to Mq. The calculation is more involved because Mq enters the 
finite parts of the renormalization constants. The calculation of P^^\g) is trivial however 
and reproduces Eq. (|9TD , i.e., P^^\g) = —g^/2Tx. Note that it also follows immediately from 
the explicit relation between the renormalized coupling and the scale Mq, Eq. ([571). Using 
Eqs. (|7l|), (|7|) and (^ we get for p^^\g): 



3a;oF^'(xo) + 2x^F^"(xo) + 2'^^^ + (Inxo + 3) 



xo- 



dxo 



(93) 



where we used 



d 



2xn 



d 



(94) 



which follows from 



dg dxo 

The calculation of 7'-^^ can be done similarly, using Z' = 1 + Z'^^y N: 



1 dz'(i) dz'(i) 



(95) 



2dlnMo dlnxo' 

To complete these calculations, we now need to specify Z'^^\ This is done with the help 
of the first of the renormalization conditions (|4^) , which has not been used so far. Using 
Eq. ( [TiD for the first derivative of V^^\ together with the expression (|77D of Z^^\ we get: 



(96) 
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Eq. (p3D becomes then: 



P'^'Kg) = -|^, f^^'Hg) = |^(lna;o + 1) (xoF^'(a;o) + x^F^"(a;o)) , (97) 

where we have written also the leading order contribution for future reference. In the small g 
limit, xq ~ oo, and from the results quoted in App. B, we get (lnxo + l)(xo-Fo (xo) + 

Xo-Fq"(xo)) ~ 1 + g'^/2n. Thus, the /3-function (|97D coincides at small g with that obtained 
earlier, Eq. (^Tj), in agreement with the usual expectation that the first two terms in the 
expansion of the /3-function are independent of the renormalization scheme provided the 
mapping between the coupling constants corresponding to the various schemes is analytic 
0]. The /5-function given by (pT]) is negative at small g, reflecting asymptotic freedom; this 
property is not affected by the corrections: (Inxo + 1) {xqFq^xq) + xIFq"{xo)) < 2, so 
j3{g) remains negative as soon as > 2. 

Finally, the function 7(1) is calculated by using Eqs. (P5|) and 



//7'(i) 

i^'\9) = xo-^ = - (xof;;{xo) + xIf;;'{xo)) . (98) 



dxo 

In the small g limit, XqFq^xo) + XqFq"(xo) ~ 1/ lna;o ~ g^/{27i). Thus, when g is small this 
expression of 7'-^^ coincides with that of Eq. (|92|) . 

One may now verify explicitly that, as claimed in the previous subsection, the expressions 
(^) of the fermion mass Mf and (^) of xm are renormalization group invariant. Consider 
first Mf. In this case, the verification can be done by a direct calculation, but it is more 
instructive to proceed by considering first the variation of Mmin = ^mi!i + (V^)^min) with 
M^°i given by Eq. (^) and M^Vn given by Eq. One finds: 

M - ^ m(°) X ^^'^'^ r99^ 



Next, from Eq. (|80|) we get: 



Note that we can replace M^j^^ by Mj in Eq. (|99D . It is then obvious that the variation of 
Mmin exactly compensates that of M^, showing that Mf = Mmin + M-^ is indeed independent 

of Mq. 

The relation (|99D can be viewed, quite generally, as a direct consequence of the invariance 
of the effective potential under renormalization group transformations. Indeed, consider such 
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a transformation, in which Mq ^ M'q, M ^ M', g ^ g' and V{M'; M'^.g') = V{M; Mo,g). 
The relation between M and M' is determined, for an infinitesimal variation dMo, by the 
renormalization group equation (^71) . The same transformation relates the values of M at 
the minimum before and after the transformation, that is. 

Mo dM^in , - , . 

= -71^) (101) 



Mmin dMo 

which, when expanded to order is Eq. (p9|) (given the relation ( p5|) for 7). 



Consider now xm- From the definition (|86D one gets 



dxM M2 1 M2 dZ'(i) , , 

-27727 + T7i722^o-^, (102) 



dlnMo M]' N Mj " rfxo 



which, according to Eq. (95) vanishes at order This may also be seen, perhaps more 

directly, by writing the relations between Mb and M as follows: 



Mb = VZ'M 



1 




A 




In In 




2N 




Wf 


1 




A 




In In 




2N 







2N 



(103) 



This equation shows that, up to terms of order the quantity (^1 + Z'^^'^ /2N^ M = 

Mfy/xM remains constant when Mo varies, with Mb and A kept fixed. 

Before closing this subsection, let us discuss the relation between the renormalized cou- 
pling constant g and the scale Mo. At leading order, this is given by Eq. (|57D. To go beyond 
leading order, we use the relation dg = f3{g)dlnMQ to write: 

The result of the numerical integration is plotted in Fig. |^, for the choice gi = ^/^T. For this 
value of gi, Mo{gi) = Mf in leading order; at order using Eq. (^21) , we get 

Moigi) = Mj (1 - (F^(l) - F^'(l) + 0) ^ M^(l - ^). (105) 

An approximate analytical evaluation can also be obtained by expanding the integrand in 

my 

1 r f, 1 f3^'^\ ..... 
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FIG. 7: The quantity ln(Mo/M/) as a function of Tr/g^ - 1 = (lnxo)/2, for N=3. Full line: 
numerical integration; dashed line: leading order relation (from Eq. (^7|)); long dashed line: 
approximate integration according to equation ( |107| ). 

By using Eq. (plf), and Xq^Qi) = 1, one then gets: 

Moig) = Mf exp |1 Inxo - (F^(xo) - a;o(lnxo + 1)^0 (^o) + O} • (107) 

This expression is identical to Eq. (^) that has been used to eliminate the scale dependence 
in the effective potential. In contrast to the linearized version, Eq. (P^, it holds even in cases 
where Mq ^ Mj, i.e., at weak coupling, where terms involving ln(Mo/Mj) may become of 
order N] such large logarithms do not enter the /3-function (as we have seen earlier, (3^^^ 
remains small at weak coupling), but can be generated by the integration over a sufficiently 
large range of values of g in Eq. (|106|) . 

The effective potential that we have obtained does not depend on the scale Mq (at order 
1/A^) , but the separation of the mass between a contribution coming from the minimum 
of the potential M^m ~ ('^V') and one coming from the self-energy S depends on the scale 
Mq. In particular, since (see Eqs. ( pSf ) and (p^)) 
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M^^/Mf may become of order unity if g is too small (see Eq. (|80|)); this, however, occurs 
only for very small such that g^ ^ e~^. 

Finally let us mention that while the choice of the scale for the effective potential is not 
an issue at this point (since the effective potential does not depend on g), we shall see that 
at high temperature the coupling constant effectively reappears in the calculation, at a scale 
determined by the temperature. This will be discussed in subsection |V D[ 



IV. THE RENORMALIZED EFFECTIVE POTENTIAL AT FINITE TEMPERA- 
TURE 

We turn now to the computation of the effective potential at finite temperature. At lead- 
ing order in the expansion, the effective potential can be split into a zero temperature 
contribution and a finite temperature one which is free of ultraviolet divergences. Thus, at 
leading order, the technical difficulties associated with the elimination of divergences are 
localized in the contribution which has been studied in the previous section. Beyond leading 
order however, the situation becomes more complicated and, at some stage of the calculation, 
ultraviolet divergences appear in terms which depend on the temperature. As mentioned 
in the introduction, one expects on general grounds 0, ||, ^ that such contributions should 
eventually cancel without the need for counterterms other than those introduced at zero 
temperature. We shall verify explicitly in this section that this indeed happens. 



A. The leading order contribution 

At leading order in the expansion, the effective potential is given by the fermionic 
contribution, i.e., the first two terms in Eq. (p5l). As we just mentioned, it can be written 
as the sum of a "zero temperature" contribution, V"*^*'-*(M), and a finite temperature one, 
\/(o)(M, T), which vanishes as T — 0. We shall denote the complete bare potential at 
leading order by V^\Mb,T), i.e., V§\Mb,T) = V^°\Mb) + V^°\Mb,T). To isolate the 



zero temperature contribution, we rewrite the sum over Matsubara frequencies in Eq. (25) 
as 

-Tyin l+ ^ ^ \ = i nculnl + 109 
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where 

n{ij) = . ^ ^ (110) 

is the fermionic statistical factor, and the integration contour is a set of circles surrounding 
each of the Matsubara frequencies. By deforming this contour into two lines parallel to 
the imaginary axis, i.e., by writing / = J^-^^^ + ji^^^^ , and using the property n{—uj) = 
1 — n{u!), one can rewrite the integral as 



-too+e 



2 



loo+e 



The term which does not contain the statistical factor gives the zero temperature contri- 
bution and it has been dealt with in subsection [III A| . The term with the statistical factor 
is the finite temperature contribution V^\Mb,T). It can be calculated by performing a 
further deformation of the integration contour so as to pick up the singularities on the real 
u-axis. We get: 

vPiMs,T) = -^ r^ln(l + e-^M. (112) 



where Ep = _^ j^^2 _ 

At this order, the renormalized potential is immediately obtained by substituting M for 
Mb in the equation above. By combining with the zero temperature contribution (|55|) , one 
obtains then the leading order renormalized effective potential at finite temperature: 

A noteworthy feature of this expression is that it is analytic around M = 0, in contrast to 
the zero temperature contribution. To see that, consider the derivative of V^'^\M,T) with 
respect to M^: 

dV^^\M,T) 1, M2 /-oo dp 1 , 



(9M2 47r M] Jo nEpel^^p + l' 

The integral in Eq. ( |114| ) is infrared divergent when M — ^ 0. By isolating the (logarithmic) 
singularity one gets [08] : 
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FIG. 8: The leading order renormalized potential V^^^ (divided by A^M|/47r) as a function of xm, 
for various temperatures T (in units of the fermion mass Mf). The temperature T = 0.57 Mf is 
the temperature Tc discussed in the next section (see Eq. ( |134| )), and 0.62Mj = l.lTc. 



where je is Euler's constant. This formula shows that the logarithms ~ InM^, at the origin 
of the non-analyticity, cancel in Eqs. ( |1 14|) and (|113|) . 

As we did in subsection [111 A| , it is convenient to express the effective potential in terms 
of the variable xm = M'^/Mj (see Eq. (p5D). With a slight abuse in the notation, one can 
then write: 



IVi X POO 

V^'^\xm,T) = -^ XM(lnxAf-l)-- / duln(l + e- 
4tt I Jo 



by/ U^+Xf. 



(116) 



where b = Mf/T. 

A plot of the leading order effective potential is given in Fig. ||. As seen on this figure, as 
T increases, the minimum is shifted to lower and lower values of xm and eventually reaches 
Xm = 0. We come back to this in subsection [V A| when discussing the restoration of chiral 
symmetry. 
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B. The 1/N contribution and the total effective potential 



As was the case already at T=0 (see subsection [1IIB|) , at T 7^ both the fermionic 
and the bosonic parts in (|25|) contribute at next-leading-order. The fermionic contribution 
originates in the next-to- leading order renormalization of Mb and Qb ™- yB\^Bi T). The 
renormalization of Vb\Mb,T) in Eq. ( |112|) amounts simply to the replacement — >■ 
M2(1 + Z'(i)/iV), with given by Eq. (0). The resulting expression can be written as 
yW(M,T) + {l/N)vj^\M,T), with V°{0){M,T) given by Eq. (|Tg) with Mb ^ M and 



where Ek = yW+M^ and = n{uj = Ej.) (see Eq. (|110|) ). The first term in Eq. (|117|) 



introduces a new divergence, not present at T=0, and which depends explicitly on the 
temperature. From the discussion at the beginning of this section, such a divergence should 
not remain in the final result, and indeed we shall see shortly that it is canceled by a similar 
one coming from the finite temperature contribution of the bosonic part, to which we now 
turn. 

The bosonic contribution comes from the last term of Eq. (p5|): 

H;1](Mb, T) = ^J [d^Q] In [1 + glUiu, q; Mb, T, A)] . (118) 

It is convenient to write the sum over Matsubara frequencies in Eq. (|118| ) as a contour 
integral, as we did for the fermionic contribution in Eq. (|109|) . By exploiting the parity 
property n(— cu, q) = Il{uj, q), one can then write: 

V,%{Mb, T) = \j^ f^2. ^ (1 + 2N{uj)) In [l + glU{u, q- Mb, T, A)] , (119) 

where 

^(^) = (120) 

is the bosonic statistical factor. 

The expression (|119|) of the bosonic part of the effective potential is analogous to the 
corresponding one for the fermionic part, Eq . ( [L 1 1| ) . However in contrast to what happened 
in the fermionic case where Eq. ( P-ll| ) leads naturally to a separation between zero tempera- 
ture and non-zero temperature contributions, this is not quite so here: the term containing 
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FIG. 9: The next-to- leading order renormalized potential V^^^ (divided by Mj/Air) as a function 
of xm > XMtiT), for various temperatures T (in units of the fermion mass Mf). The temperature 
T = 0.57Mf is the temperature Tc discussed in the next section (see Eq Jl34D ), and 0.62Mj = l.lTc. 



the statistical factor N{uj) certainly vanishes as T — > 0, and as such can be considered as 
a genuine finite temperature contribution; the term without the statistical factor reduces 
to the zero temperature effective potential as T 0, but it does contain finite temper- 
ature contributions (since 11 depends now on T). Nevertheless, the separation suggested 
by Eg . (11191) is useful as it allows us to clearly isolate the ultraviolet divergences: these are 
contained entirely in the term without the statistical factor. We shall then add the zero 
temperature counterterms (|^) to that term, renormalize the mass and the coupling con- 
stant with the leading order renormalization constants (Z^^^ from Eq. ( ^1]) and Z'^^^ = 1), 
and rewrite the renormalized effective potential V^^'^ (M, T) as the sum of two contributions: 
vi'\M,T) = vi]\M,T) + V,^'\m,T) with 



vff{M,T) 



1 /-A d\ 



(2vr) 



In 



D-\zqo,q;M,T) 
D,\QE;Mf) 



+ 



M2 A2 
^^In — ^ H 

Mj 4:71 



4tt 



M2 a 
— - In In — 

M 



121) 



/ 



WJ(M,T) 



dq 



N{ij) In \D-\u;,q;M,T) 



122) 



-oo 27r J — ioo+e 27rz 

where D~^{iqQ,q; M,T) = Dq^{Qe',M) + g'^Il{iqQ,q; M,T) is the renormalized inverse a- 
propagator studied in App. A. 
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Clearly, when T ^ 0, vl;^^{M,T) becomes the zero temperature part of the bosonic 
contribution to the potential analyzed in subsection |111 B| ; in the same limit 14^2^ (M,T) 



vanishes. The divergent contributions to Vj,''^^(M, T = 0) are dealt with as in subsection |111 B| . 
When the temperature is non zero, a new divergence occurs in v/}''(M, T). This can be 
isolated with the help of the asymptotic behaviors of Dq^ and II given in App. A (Eqs. 
and ( [A26| ), respectively). One gets: 

vS(M,T)~-MH.M^)ir^4'^). (123) 



This divergence cancels that of the renormalized fermionic contribution vi^^ of Eq. ( |117D , 



as anticipated (the two A-dependent terms in Eqs. (117) and ( |123[ ) are identical to within 
terms which vanish as A — 00). 

At this point, it is useful to go one step further and, in analogy with what we did in 
Eq. (^) of the previous section, we absorb the finite part of the counterterm into the 
redefinition of the variable xm according to Eq. (|86|). Adding V^^\M,T) and the second 
term in Eq. ( [117| ) one can write (discarding terms of order 1/N'^): 

/3 io TT V J N Jo -K Ei, b Jo Tc \ r ' 

with Xm = M^/M] (1 + Z'^^^/N), as defined in Eq. (|6|), and b = Mj/T. Thus, in 
complete analogy with the zero temperature case (see Eq. (P7|)), after the redefinition of the 
variable xm, the leading order contribution to the renormalized potential keeps the form of 
V^"^{xm,T) given by Eq. ([TTBI) . 



By using the results just obtained and defining new functions F{xm,T), and G{xm,T), 
we can then write the complete effective potential at order in the form: 



VixM,T) = -J- 



8 a"t 

ArxAf(lnxA/ - 1) - 7 / duln(l + e-^V" 

b Jo 

+e(xM - 1) + FixM, T) + G(XM, T)l , (125) 



where the functions F and G are given in App. B. The function F{x,T) is an extension 
at finite temperature of the function Fq{x) introduced in the previous section; to within an 
additive constant, it is proportional to the finite part of V^^i{M,T). The function G{x,T) 
is proportional to V^^ {M,T): 

=^GixM,T). (126) 
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It can be calculated by deforming the integration contour so that it goes around the 
singularities on the real axis. This operation requires that grand circles at infinity generate 
no contribution, and that the discontinuity of the integrand vanishes at = 0, properties 
that are both satisfied. The resulting integral along the cut is then finite because of the 
statistical factor. It can be written as: 



dg dcj 

■ ;v I iTi 

D-\uJ -ie,q;M,T) 

+00 dq roc 



J -00 Ztt Jo Zm 



D-\u; + ie,q;M, T) 



r + 00 (^Q fOO (^ij 

= TT- —N{uj)2lmlnD-\uJ + te,q;M,T). (127) 

j-00 2tt Jo 2tc 

One may also express the argument of the logarithm in terms of the phase shift 6 of 
quark-quark scattering (see e.g. Eq. ( pi6| )): 



D-\uj + te,q;M,T) 



~ D-^iLU~te,q;M,Ty ^ ^' 



with 



2 lmn.{uj + ie,q;M,T) 



tan S(uj, q; M, T) = ' (129) 

^ ' ^ ^ ReD-^{u} + ie,q;M,T) ^ ' 

The latter are the expressions we use to calculate the function G{xm, T) numerically and to 
understand its behavior for small values of xm (see App. B). 

To get an intuitive picture for what these formulae represent, let us first rewrite the 
expression (|127|) by integrating by parts, 

(The boundary term does not contribute because 5{uj, q) vanishes at small u in the same 
way as Imn(u;, q).) Then, let us imagine that the a excitation corresponds to a simple pole 
of the propagator, i.e., assume D''^{u, q) = uj"^ — oj"^. In this case, d5 jdu = 7r6{uj — Ug), and 
V^^^^{M) is simply 

which is the free energy density of noninterating bosons with energies Ug. In the present 
case however, the sigma excitation does not correspond to a pole of the propagator (see 
subsection |VB|) , and this simple picture does not quite hold. 



The 1/A^ contribution to the renormalized effective potential given by the second line in 
Eq. (|125|) is displayed in Fig. |. As was the case at zero temperature (see the remark at the 
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FIG. 10: The effective potential (divided by NMf/iir) as a function of xm > x m, {T) , for various 
values of N and the temperature T. The temperature T = 0.57Mf is the temperature Tc discussed 
in the next section (see Eg. [134]) ), and 0.62Mj = l.lTc. 



end of subsection B| )), there is a minimum value xa/, (T) of xm below which a Landau 
pole appears and the function F{xm,T) becomes complex. For this reason, the next-to- 
leading order contribution to the potential is plotted in Fig. |] only for xm > xm^(T) (the 
way Xm, depends on the temperature is discussed in subsection [VB|) . In the same region 
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xm < XM^iT) where F is not real, the function G takes anomalously large and negative 
values (see Fig. ^ and the discussion in App. B). For xm ~ xm,, both G{xm,T) and 
F{xm,T), and therefore also Vi{xm), have a large slope as can be seen in Fig. ||. 

The complete effective potential is displayed in Fig. |Ty for various values of and the 
temperature T, as a function of xm (for xm > xm,{T))- The curves labelled = 100 are 
very close to those corresponding to the leading order potential (not drawn). If we follow 
these curves, one recovers the behavior of the leading order potential, already discussed: 
as the temperature increases, the minimum of the potential shifts to lower values of xm 
and becomes more and more shallow, until it reaches xm = where it becomes again a 



pronounced minimum. A striking feature of the curves in Fig. |T0| is the large effect of the 
1/A^ corrections: already at moderate temperatures, one sees that the minimum disappears 
for small values of A^. This is to be related to the large slopes at small xm of the curves 
in Fig. which drive quickly the minimum to low values of Xm-, making this minimum 
dangerously close to the threshold xm^ for the appearance of the Landau pole. Note however 
that, at larger temperatures, the potential has a well defined minimum at xm = for all 
values of A^. 

All these properties of the effective potential, and their physical interpretation, will be 
discussed at length in the next section. 

V. THERMODYNAMICS 

We use now the results that have been established in the previous sections to discuss the 
thermodynamical properties of the system. We start, in the next subsection, by reviewing 
known leading order results concerning the restoration of chiral symmetry at a finite tem- 
perature Tc. In this one-dimensional system, this phenomenon is specific of the mean field, 
or large A^, approximation based on uniform quark condensates, i.e., on constant solutions 
of the gap equation (p^f). (As argued in Ref. |T^, kink configurations of the condensate, if 
taken into account, would presumably prevent chiral symmetry breaking at any non-zero 
temperature). At next-to-leading order in the 1/A^ expansion the fluctuations around the 
uniform condensate provide a small correction to the mean field picture as long as the tem- 
perature remains small. But as the temperature increases these fluctuations become large. 



eventually making the behavior of the system pathological; as we shall see in subsection V_C, 
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this results in a breakdown of the expansion at some temperature below the mean field 
transition temperature Tc. We shall see that this breakdown is related to the presence of a 
Landau pole which is discussed in subsection |V B. 



The last two subsections are devoted to the high temperature limit (T ^ T^). There 
chiral symmetry is realized, and the expansion remains a good approximation scheme. 
It provides a description of the system as weakly interacting fermions, as expected from 
asymptotic freedom. Quite remarkably, while the effective potential obtained in the previous 
section does not depend explicitly on the coupling constant, we shall see that the pressure 
can nevertheless be expanded in terms of an effective coupling whose magnitude decreases 
logarithmically with the temperature. To better understand this result, a direct perturbative 
calculation of the pressure is presented in the last subsection. 



A. Restoration of chiral symmetry in leading order 

The finite temperature effective potential at leading order may be written (see Eq. (|116|) ): 



with Ek = vP + M^. As discussed in subsection IIC, at this order the minimum of Eq. 
( |132| ) is simply the quark mass Mf{T). It is a decreasing function of T. A simple analysis 
shows that, at small T, this decrease is exponentially small: 



e~''f/'. (133) 



As the temperature increases further, however, Mf{T) becomes smaller and smaller and 
eventually vanishes, as can be seen from the plot in Fig. |lT] which displays the function 
Mf{T) obtained numerically. One may use the expansion in Eq. (|1 1 5|) to calculate the 
temperature Tc at which Mf{T) vanishes |]14|, ITSf. One gets: 



Tc = Mf ^ 0.567M/. (134) 

vr 

As T —> Tc the system undergoes a second order phase transition with mean field expo- 
nents. To verify this explicitly, we expand the effective potential around M = 0: 

V^(o) (M, T) ^ V W (0, T) + a + 6 ' (135) 
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FIG. 11: The fermion mass Mf{T) in leading order in the 1/A^ expansion as a function of T/Mf. 
Also indicated is the value M^:(T) of M below which a Landau pole appears (see subsection |VB| 
and App. A). 

with 

M| ^^TCmfl (136) 

where C(3) ~ 1.202. The coefficient a follows easily form the expansion ( |115|) . The coefficient 
h is obtained by using the expression ( p.09| ) of the potential as a sum over the Matsubara 
frequencies. By taking the second derivative of this expression with respect to one 
obtains a convergent sum which is proportional to h: 

1 _ 7C(3) _ 2b 

A ^ [\2n + l\7cTf 167r3T2 Mj' ^ ' 

(Note that the expansion ( p.35| ) of the potential is meaningful only at finite temperature, as 
is obvious on the form of the coefficients a and h which are singular in the limit T 0.) 



The form (|135|) of the effective potential controls the behavior of Mf{T) / Mf in the vicinity 



of Tc- A simple calculation gives: 



MfiT) -a 



^' /^-^ (138) 



Mf V 26 \ 7C(3)M/V 
Above Tc, the fermion mass Mf{T) vanishes, and so does the (renormalized) quark con- 
densate: ('?/''?/') = —NMf{T)/g'^ (see Eq. (p9|) ). Thus, at Tc the discrete chiral symmetry is 
restored. 
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The pressure is simply related to the effective potential at its minimum: P(T) = NP^^^ = 
—NV^^\M = Mf{T),T). The entropy density can be easily calculated from s{T) = dP/dT. 
Since V^^\M, T) is stationary for variations of M around the value M = Mf(T), the deriva- 
tive can be simply taken at fixed M. The result is that the entropy density is just that of 
free particles of mass Mf{T): 

/oo fjU 
— [(1 - rik) ln(l - rik) - In n^] , (139) 
-oo ZTT 



where = n{Ek) is the fermion statistical factor (see Eq. (|110|) ), and = '^k? + M'j{T). 
When T > T^, the fermion mass vanishes and we have (/q°° dx\n{l + e~^) = 7r^/12): 

p,o, ^ ^ ,,0, ^iP^^^^S. (140) 

6 dT 3 

Thus for T > T^, interaction effects completely disappear, and the system behaves as a 
system of free massless fermions. 

We show in the last two subsections that the corrections do not alter much this 
picture at high temperature where, due to asymptotic freedom, the effect of the interactions 
remains small. But the corrections do affect very strongly the behavior of the system 
below Tc, as we shall see shortly. Before turning to these corrections, we shall, in 
the next subsection, continue to discuss physical properties of the system in the mean field 
approximation; some of these will be useful in the interpretation of the next-to-leading order 
results. 



B. The a excitation and quark-quark scattering amplitude 

At zero temperature, the propagator Do{Q; Mf), obtained by setting M = Mf in 
Eq. (|X6|), describes cr-meson excitations with mass M„ = 2Mf. This is easily seen by 
solving the equation Dq^{Q^ = — M^; Mj) = for M„. Note however that this point where 
Dq^ vanishes does not correspond to a simple pole in the a propagator, but to a branch point 
limiting the region of phase space where the a excitation can decay into quark- ant iquark 
pairs. 

In order to see that we need to perform an analytic continuation from the Euclidean 
momentum Qe = {qo,q) to the Minkowski one Qm = {^^, q), with go ~iui. From Eq. ([A 61), 
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FIG. 12: The real part of D^^{uj,q;M) as a function of {uj^ - q^)/M^ for M < Mf. For M = 
Mf, RcDq^ vanishes at the single point (w^ — q^)/M'j = 4, beyond which the imaginary part 
starts to grow. For an arbitrary M, ReD(^"'^(a;, (7; M) has a minimum at this point, with value 
(^V^) ln(M/M/) (see App. A). When M = = M//e, Re D^^ vanishes at = 



one gets: 



D^\u,q-M) = ^ 

ATI 



(141) 



Let us set ujq = y + AMj. The inverse propagator Dq^{uj, q; Mf) has cuts on the real u 
axis for u > Ug and u < —ujq] the corresponding imaginary part is (for > 0): 

-,2 



(142) 



2 ^ - g2 

The point = cjg, is the branch point corresponding to the threshold mentioned above. 
A plot of ReD^^{u + ze, g; M) = Rel^o ^(o; + ie, g; Mj) + (fi-V^i") ln(M/M/) is displayed in 
Fig. (0). 

At finite temperature the threshold remains located at 2Mj(T), as long as T < Tc. To 
see that, consider the equation which defines M„] 

D-\uj = M,,q = 0;Mj{T)) = 0, (143) 

where = Dq^ + g^Il is the finite temperature a propagator studied in App. A. Let us 
set Mo- = 2Mf{T) and verify that indeed this is the solution. Using the explicit expressions 
for Il{uj, q) given in App. A, we get: 



D-\u = M„,q = 0; Mf{T)) = ^ In 



'Mf(T)\ . f 



TC 



°° dk Uk 
vr Ek 



(144) 
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with Uk as in Eq. ( p.39| ) above. This equation is precisely that which determines Mf{T), 
i.e., the gap equation obtained by equating the right hand side of Eq. ( |114| ) to zero, which 
proves the point. 

A more complete description of the a excitation is obtained from the corresponding 
spectral function p{u!,q). At zero temperature, we define {u real): 

Po{uj, q) = ^[Do{lj + ie, q; Mf) - Dq{uj - le, q; Mj)] = 21mDo{uJ + ie, q; Mf). (145) 

This is an odd function of uj. Using the formulae given above, one easily obtains (for > 0): 

4 



9 \ 



* 1 + 4 Arctanh 



At fixed g, this is a decreasing function of cj, peaked at ~ ujq\ when uj — > oo, pQ{uj,q) 
vanishes as pQ{u!,q) ~ Att"^ / {g In u"^)"^ . By replacing, in Eq. ( |145| ), Dq by the finite temper- 
ature propagator D{uj,q; Mf(T)), one obtains the spectral function at finite temperature 
p{uj,q]T). The typical behavior of p{uj,q;T) for temperatures below and above Tc is dis- 



played in Fig. |T3[ We note, below Tc, the threshold at = yg^ + 4Mj(T)2, beyond which 
the spectral function resembles po{uj,q) given by Eq. (|146|) . For < < g, a typical finite 
temperature contribution appears, that of the scattering of quark-antiquark pairs on quarks 
or antiquarks present in the heat bath (see App. A for more details). As T approaches Tc 
from below, Mf{T) — > 0, and the two contributions merge, leaving a sharp peak at small 
u. The peak structure survives just above (see the curve labelled T = 0.6 in Fig. O), 



but as one raises the temperature the spectral weight becomes spread over a larger and 
larger energy range, as can be seen in Fig. |T^. The behavior of the spectral function is then 
governed by the high temperature regime of the fermion loop discussed in App. A. In this 
regime and for large values of g (g ^ Tc), the spectral function becomes a step function at 
UJ = q followed by a long logarithmic tail. 

Consider finally the quark-quark scattering amplitude T(c(j, g; M, T), where the mass M is 
treated as an independent parameter. This is the generalization at finite temperature of the 
quantity introduced in Eq. (|6^) , i.e.: T{uj,q; M,T) = g'^D{uj,q; M,T). At zero momentum, 
we have (see Eqs. (|27|)): 



'd2y(o)(M,T)' 



-1 



T(0,0;M,T) = (7'limD(0,g;M)= dM^ J ' ^^^^^ 
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FIG. 13: The spectral function of the a excitation above Tc (top, q = 0.01) and below Tc (bottom, 
q = 0.1), as a function of co expressed in units of Mj (upper curves) or Mf{T) lower curves. 



This relation may be verified by calculating the second derivative from the expression (|113|) 
of the potential and using the results given in App. A to obtain: 



lim 



M 



- 1 + ln— +2/ — 

TT \ Mf I Jo TT 



^ dk / dn 



K n 

El dE, + 



(148) 



9-0 r (0, g; M, T) vi y ivif ^ 

The second term in the r.h.s. of Eq. ( |148| ) is limg^o n(0, g; M, T) (where we have used 
Eq. ( |A9D ) . The two terms in the integral correspond respectively to the contributions of 
scattering processes and pair creations (see App. A). As mentioned already, at finite tem- 
perature, and when M 7^ 0, the double limit ^ 0, g ^ of D{ijj, q) is not regular and the 
result depends on the order in which the two limits are taken. If the limit g — > had been 
taken first, the scattering term, proportional to dn/dE, would be absent in Eq. ( |148|) , as 
may be seen from Eq. (|A9|) . 

When M — i> 0, the scattering term vanishes and the limit — > 0, g of T{uj, g) 
becomes regular. In that same limit where M — * 0, the second term in the integral of 
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Eq. ( |148| ) (i.e., the pair creation term) develops a logarithmic divergence, which in fact 
cancels the logarithm of the T = contribution (the first term of Eq. ( |148| ) ) . Thus one has 
(see Eq. PT9|) ): 

lim r(0, 0; M, T) = , (149) 
Af^o ^ ' ' ' ^ ln(T/Tc) ^ ^ 

which is nothing but the square of the renormalized coupling constant evaluated at a scale 
Mo ~ T (see Eq. (|57|) ). Thus at high temperature, i.e., when T 3> Tc, the scattering 
amplitude becomes small, as expected from asymptotic freedom. It can be also verified that 
no Landau pole occurs in this regime. 

Below Tc the situation is complicated by the presence of a Landau pole in the a propa- 
gator. The threshold M^{T) for the appearance of the Landau pole, which coincides with 
the value of M at which the second derivative of V^'^'> vanishes, is plotted as a function of 



the temperature in Fig. |Tl|. As the temperature increases, M^, is shifted first to slightly 
larger values: this is because at moderate temperatures the finite temperature contribution 
in Eq. (148) is dominated by the first, negative, term in the integral. For temperatures 



beyond approximately Tc/2 however the second, positive, contribution becomes dominant, 
resulting in a shift of M* to lower and lower values, reaching M^, = at Tc. That M^, vanishes 
precisely at Tc follows from the simple fact that is located between two extrema of the 
potential, which merge at T^. Thus, as T approaches Tc, the fermion mass Mf{T) becomes 
close to Mj,{T), and the calculation of the next-to-leading contribution to the effective po- 
tential becomes ill defined in the vicinity of its minimum: in particular F{xm',T) becomes 
complex for xm < xm, {^m, = (^*/^/)^) and both F{xm',T) and G{xm',T) present large 
negative slopes for xm ^ xm, (see discussion in subsection [IV B| and App. B). The latter 
has important consequences which will become clear in the next subsection. 

C. The quark condensate at order 1/N 

In this subsection, we examine how the physics of chiral symmetry breaking and its 
restoration, which we have just discussed, is modified by the 1/A^ corrections. Chiral sym- 
metry breaking is characterized by a non-vanishing value of the quark-condensate, and 
the latter is proportional to the value Mmin of the minimum of the effective potential (see 



Eq. (poD). The effective potential has been obtained in subsection |IV B| to order 1/A^ in terms 
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of the renormalization group invariant xm defined in Eq. (0). We shall analyze here the 
variation with the temperature of Xmin = ^Mmin o*^ obtained at next-to-leading-order. 

The quantity XmmiT) may be determined to accuracy by following the strategy 
exposed in subsection y B. Generalizing Eq. ( [55| ) at finite temperature, we set 

where a;J^L(^) is the minimum of V^^\xm,T) in Eq. (|116D. Clearly, x^^\^{T)M] = M]{T), 
where Mf{T) is the temperature dependent quark mass introduced in subsection |V A| , and 
plotted in Fig. |lT[ The correction x[^]„(T) is given by (see Eq. (pH])): 



dy(l)(a:M,T) 




dxM 


(T) 






dxlj 





(151) 



with 

dVW(xM,T) _ ^ ^ F'rfL(r),T) + G'(xSL(r), T), (152) 



dXAf 

and d2l/(°)(xM, T)/dx|^ can be determined from the results of subsection [V A| . In Eq. (|152D 



F' and G" are the derivatives of the functions F and G with respect to xm\ they are calculated 
numerically. 



A plot of Xmin{T) is giveu in Fig. |T^. At small T, Xmm > 1, in agreement with the study 
of the effective potential at zero temperature (see the caption of Fig. ^). However, when 
the temperature T increases, the 1/A^ correction eventually turns negative at T ~ Tc/2 and 
becomes large when the temperature increases further. Thus, for instance, la^mL/^mLl — 0-1 
when T/Tc ~ 0.55, |a;^]„/a;^]„| ~ 1 when T/Tc ~ 0.65, |a;^]„/a;^]„| ~ 10 when T/Tc ~ 0.85 
and |a;^L/3^mLl — 30 when T/Tc ~ 0.92. Thus, when T/Tc ~ 0.65 the 1/A^ correction 
becomes too large to be really considered as a "correction", and the 1/A^ expansion breaks 



down. This is to be contrasted with the situation at zero temperature where 



(1) /JO) 



mm I 



0.25 (see Eq. (|88D and the discussion after). Note that what makes la^mln/^mLl large as T 
approaches Tc is the numerator dV^^ydxM (the denominator d'^V^'^y dxl,^ goes to a constant 
as T — i> Tc); as we have already discussed, both F\xm) and G'{xm) become large and 
negative as xm approaches xm, , and Xmin does approach xm, as the temperature increases. 

This pathological result is confirmed by a direct (numerical) minimization of the potential 
in Eq. ( |125| ). For all values of and T, one finds values of Xmin{T) very close to those 
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1.2 n 




T/Mf 

FIG. 14: The value of 

x-min calculated from Eqs. ( |15[1| ) and ( |151| ), as a function of the temperature 
T (in units of Mf), and for various values of N. The curve labelled L.O. corresponds to the leading 
order potential. 

deduced from Eqs. (|150|) - (|151|) whenever the minimum exists (the two evaluations of Xmin{T) 
may differ a priori by terms of order l/N"^). However, as anticipated in subsection |IVB| , the 
potential no longer has a minimum when T exceeds a certain value, which depends on A^: 
for = 100, the minimum ceases to exist when T > 0.93Tc, for = 10 when T > 0.63Tc 
and for = 3 when T > O.SIT^. Again the disappearance of the minimum may be related 
to the rapid drop of the functions F{xm,T) and G{xm,T) in the vicinity of xm*, the value 
of xm below which a Landau pole appears, as discussed in App. B. 

At this point one could speculate and try to relate the breakdown of the 1/A^ expansion 
to the fact that the true transition temperature of the model is presumably Tc = ||TU| . The 
mean field approximation which ignores part of the important degrees of freedom (e.g. the 
kink configurations) gives a poor representation of the physics of the system at finite tem- 
perature. This mean field physics persists for moderate temperatures: then the fluctuations 
around the uniform condensate do not change significantly the state of the system. But be- 
yond some value of the temperature the fluctuations become dominant and perhaps mimic 
the effect of degrees of freedom left-out in the mean field approximation: these fluctuations 
are responsible for the rapid decrease of Xmin for T ^ Tc/2 which drives the system towards 
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T/Mf 



FIG. 15: The pressure as a function of T/Mj. The full line is the pressure of free massless fermions, 
Eq. ( |14[1| ). The thick dashed line is the pressure at leading order (L.O.). The other curves include 
the corrections, for = 10, = 5 and = 3. For small temperatures, the pressure is an 
increasing function of A^; the curves corresponding to the various values of A^ cross each other for 
T ^ 0.5Mj, and for larger temperature, the pressure decreases with increasing A^. 

its true equilibrium state, where chiral symmetry is restored. 

Another signal of the inadequacy of the 1/A^ expansion for temperatures T ^ Tc is 



provided by the study of the pressure. A plot of this quantity is given in Fig. Note that 
although the determination of the minimum of the effective potential becomes inaccurate 
in the vicinity of Tc (and may not even exist), this is of little consequence for the present 
calculation: since the potential is flat in the vicinity of Tc, its value, i.e. the pressure, is 
well estimated. One sees in particular that the entropy density, i.e., the derivative of the 
pressure with respect to T, smoothly increases with the temperature below Tc, but decreases 
suddenly at Tc, suggesting a first order transition. This is clearly unphysical. In fact the 
plot reveals two regimes: at low temperature the system exhibits unphysical mean field 
behavior; this stops abruptly around Tc where a new regime sets in. This regime is that 
of high temperature with restored chiral symmetry: there, the approximation appears 



51 



to be a good approximation; it allows us to study the physics of the system in conditions 
where, because of asymptotic freedom, it is expected to behave as a gas of weakly interacting 
massless fermions. We shall verify in the last two subsections that this is indeed true. 



D. Thermodynamics at high temperature 

We discuss now the limit of high temperature, T ^ T^, where the condensate vanishes 
and, in leading order, the quarks are massless. As we shall see, the thermodynamical 
functions can then be expanded in powers of 7r/ln(T/Tc) which we shall interpret as the 
running coupling constant at a scale of the order of the temperature. This provides a nice, 
and non trivial, illustration of the behavior expected from asymptotic freedom. Let us 
consider the limit of the contribution to the effective potential V^-^\ i.e., the second 
line in Eq. ( |125| ) when T ^ T^. It is obtained from the high temperature limit of F{0, T) + 
G{0, T), the constant C, playing no role in this limit. The function F{0, T) is given by the first 
line of Eq. (|B15|) (the second line is a finite constant and plays no role at high temperature). 
It can be written as 



F(0,T) = —2- / ds dr In 



Mj 71 Jo 







1 + g In (T2/M| ^{s^ + r2)2 + Mf/T^) 
where the quantity g that we have introduced in this equation is 



(153) 



^2 ^ 



which coincides with the scattering amplitude at zero momentum for massless particles (see 
Eq. ( p.49| )). According to Eq. (|57D, this quantity may also be interpreted as the effective 
coupling constant at the scale Mq ~ T, i.e., g"^ = g'^{MQ = t^T / e""^^^). We shall come back 
to this identification in the next subsection. As is clear on Eq. (|153|), aside from the overall 



factor in front of the integral, all the temperature dependence of the effective potential is 
entirely contained in g^ (the extra temperature dependence in the denominator in Eq. (|153|) 
is numerically negligible). 

Since g becomes small as T increases, one may attempt to expand V^*'P(T) in powers 
of g. In order to do so, we note that Bxiis, r) — In (^{s^ + r'^)T^/Mj^ decreases rapidly as 
r,s ^ 1 (see App. A), and the contribution of this region to the integrals is negligeable. 
The integrand is also regular in the r, s — limit, and correspondingly the contribution to 
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integrals of the region r, s <C Mj/T"^ is negligeable. We conclude that the only region of the 
s, r plane which contributes to the integrals in the limit T ^ Tc is the region r, s ~ 1. There, 
both Bx{is, r) — In (^{s"^ + r'^)T^/Mj^ and ln(s^ + r^) are small compared with ln(T/Tc) when 
T ^ Tc. One can then proceed to an expansion of F{0, T) to get 



T2 /g2 g4 



(155) 



with 



2 

vr Jo 



ds 



dr 



BT{is,r)-\Yii{s' + r')T^/Mj 



(156) 



and 



F, = — 



TT JO 



ds dr f r)]^ - [in ((s^ + P)T^/M^f 



(157) 



F2 can be calculated analytically by noting that Bx{is,r) - ln((s2 + r2)T2/M2 
2TTll{iTs,Tr; M = 0,T), that is, F2 is proportional to the (go,g) integral of n(igo,Q';^ = 
0,T). It can be calculated, by using Eq. ( |X9| ) for 11, performing first the go integral, 
changing to the variables k± = k ± q/2 and doing the k± integrals, and finally using 
dyln{l + e~y) = 7r^/12 for the remaining integration. F4 has been calculated numer- 
ically. We obtain then: 

2 



-TT^ ^ 6.5797 Fa ^ -4.322 



(158) 



Let us now turn to the evaluation of the function G{0, T), given by Eq. ( [B16D . It can be 
written as: 



G(0,T) 



M]7f 



ds 



dr N{s)5{s,r;T) 



(159) 



where s = uj/T, r = q/T, N{s) = l/{e^ — 1), and: 

lmBT{s, r) 



tan (5 



(160) 



ln(rVT2) + ReST(s,r)' 

To obtain the limit of G(0, T) forT ^ Tc one observes that the statistical factor N{s) 
limits the s-integral to s ^ 2. For these values of s, the imaginary part is significant only 
in a finite range of values of r (e.g. r ^ 10 for s ~ 1), and vanishes exponentially at larger 
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r. In this domain of values of s and r, the factor IieBT{s,r) in the denominator remains of 
order 1 while the term ln(T/Tc) becomes large when T ^ Tc. One can then expand: 

~2 

. X j^lmBTis,r) 9^ ^ ^ 

tano = — — — r~ ^ [mBT{s,r) 

l + ^BT{s,r) 27r 

The expansion of G{0, T) follows: 
where 

/-oo /-oo 

G2 = - dr dsN{s)lmBT{s,r), (163) 

TT -/O Jo 

and 

A /-oo /-oo 

= / rfr / rfsAr(s)ImST(s,r)ReST(s,r). (164) 

71 Jo Jo 

G2 can be calculated analytically, by noting that lmBT{s,r) = 27rIm(n(Ts, Tr, T)). Using 
then Eq. ( [A22[ ) , we can perform the two integrals using dyln{l — e~^) = 7r^/6. G4 has 
been calculated numerically. We obtain: 

G2 = --TT^ ^ -6.5797, G4 ~ 6.281 . (165) 
3 

At this stage, several comments are in order. Recall that, at leading order, the pres- 
sure is that of non interacting massless particles as soon as T > T^. P^^^ = (see 
Eq. ( |14CI|) ); thus, when T > T^. all the effects of the interactions are contained in V^^^-*. 
Asymptotic freedom leads us to expect that, at high temperature, these interactions should 
be weak. Since the coupling constant has completely disappeared in the final expression of 
V'^^'> this property is not immediately obvious, but the analysis of this subsection indicates 
how this happens. The next-to-leading-order contribution to the pressure can be written as 
p(i)(7') = T^/(^^) where, aside from the explicit which carries the dimension, (almost) 
all temperature dependence is contained in the effective coupling g. This effective constant 
is nothing but the running coupling at a scale ~ T, and it becomes small at large T. The 
function /(^^), which vanishes for g'^ = 0, can be expanded in powers of g, and we have given 
above the leading contributions of this expansion. A noteworthy feature of this expansion 
is the vanishing of the term of order g"^: F2 + G2 = (see Eqs. ( p.58| ) and ( [L63| )). We shall 
clarify in the next subsection the origin of this cancellation and, more generally, verify that 
the properties of the pressure at high temperature that emerge from the 1/A^ expansion are 
those one expects form perturbation theory. 
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l_|_ReEr(s,r) + 



(161) 



E. Perturbation theory at high temperature 



We shall now reconsider the high temperature limit within ordinary perturbation theory, 
i.e. within an expansion in powers of the coupling constant starting directly from the 
lagrangian (|1]). We shall then be able to recover the results of the previous subsection, and 
understand the origin of the properties of the pressure that we have just discussed. Note 
that in contrast to what happens at zero temperature, where perturbation theory about 
the symmetric vacuum is meaningless because of infrared divergences, at finite temperature 
no such divergences occur because the fermions are effectively massive (because of their 
non vanishing Matsubara frequencies). Consider then the diagrams contributing to the 
pressure up to order g'^ (recall that the pressure is minus the effective potential). These 
are displayed in Fig. |16|. Note that these are all included in the calculation of the pressure 
at order That is, when expanding the result of our calculation in powers of g, 

one should reproduce the perturbative expansion up to order g^. Let us call Pq, P^, Pc, Pd 
their respective contributions. To handle their ultraviolet divergences, we shall follow the 
standard procedure of perturbation theory. This implies using a renormalization scheme 
somewhat different from that of the rest of this paper, but this has no consequence for the 
final results. In the first part of our analysis, we shall assume that the quarks have a small 
constant mass m, i.e. a term mipil) is added to the lagrangian (|T]). This mass will serve both 
to avoid spurious infrared divergences in intermediate calculations, and also to keep finite 
expressions which would vanish if m were zero to start with (in fact, all contributions except 
Pd vanish when m ^ 0): keeping m finite is then useful to illustrate the systematics of the 
cancellation of ultraviolet divergences and of the various renormalizations involved. Once 
this will be understood we shall take the limit m — >■ 0. 

Consider then the first contribution: 



Pa 



dk m ,^ ^ . 
- — (l-2n,) 



(166) 



where Ek = y/k'^ + m? and = n{Ek) is the fermion statistical factor (|110|) . Note that 
Pa is proportional to {ipip)"^ (given by Eq. (pIS]) with M ^ m), from which one could guess 
that Pa vanishes in the chiral limit m = 0. However, as it stands. Pa is ill defined and needs 
renormalisation. To do that, we rewrite Pa as follows: 



Pa 



dk m 
2tt^ 



dk m 
2tt^ 



dk m 
27tE 



-2ni 



+ 



dk m 

277 



(167) 
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(a) 



(b) 
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(c) 



(d) 



FIG. 16: The diagrams of perturbation theory contributing to the pressure to order ((a) and 
(b)) and ((c) and (d)). The dotted hne represents the fermionic interaction (proportional to 
g'^/N), not the fuh fi-propagator. The diagrams (a) and (c) are of order N, the diagrams (b) and 
(d) of order 1. 

The first term in the second line of Eq. ( |1671 ) is a vacuum contribution and it can be 
discarded. The last term is finite. The potential difficulty comes from the second term which 
is divergent and temperature dependent. However such a divergence can be eliminated by a 
mass renormalization, i.e., by adding to the Lagrangian the counterterm 6m2ipip with 6m2 
chosen so as to cancel the correction of the zero temperature mass at order g'^, i.e., 




(168) 



With the mass counterterm included. Pa becomes 




(169) 



and the unwanted contribution disappears in P^, as expected. 
Consider next the "exchange" term 
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One can notice that n(r = 0, x = 0; m, T) = lim^^o TrS'(r, x = 0;m, T)S{—t, x = 0;m,T), 
where S is the fermion propagator (see Eq. (|C17| )), to write Eq. ( |17(]| ) as: 

1 + kq 



;i-2nfe)(l-2n, 



(171) 



2EkEq J 

Note that before taking the hmit t 0, the integrand in Eq. ( |171D contains exponential 
factors e^'^*'«'^ which, together with statistical factors, insure the convergence of the integrals 
over k and q. In the limit r = 0, the first term in Eq. ( |1 71|) is a divergent constant 
independent of the temperature and it can be ignored. The contribution proportional to kq 
can be written as a product of two integrals which vanish by parity (at any finite r). What 
remains then in Eq. ( |1 70| ) is: 

^ '■ dA; m , . _ ,\ f p dq m , ^ _ A Pn 



-ill 



2nr, 



2N' 



(172) 



Thus the exchange term is just proportional to the direct contribution Pa calculated above. 
It can be renormalized in the same way. We note now that all the contributions considered 
so far remain, after renormalization, proportional to m. It follows that all the contributions 
of order g"^ to the pressure cancel in the limit m —>■ 0: such contributions are proportional 
to {ipip)"^, a quantity which vanishes when m = 0. 

Let us turn now to the first of the order contributions. 



g^N 



dk m 



2ni 



n(0,0; m,T) = -g^PaU{0,0;m,T) 



173) 



where n(0,0;m, T) stands for Il{uj = 0,q = 0;m, T). One divergence contained in this 
quantity is eliminated by coupling constant renormalization. The needed counterterm is 
obtained form Pa of Eq. ( |167| ) in which one replaces g'^ by g^ + 5''^no(0, 0; Mq). One gets 
then: 



P" 



dk m 

27T^ 



2ni 



:i + /no(0,0;Mo)) 



;i74) 



By combining P^ and the order g'^ of P^', one gets 



P' 



g^N 



dk m 
2^E'J 



2nk 



(n(0,0;m,T)-no(0,0;Mo)) 



;i75) 



where the quantity 11(0, 0; m, T)— no(0, 0; Mq) is finite. The contribution is still divergent, 
but the divergences can be eliminated by the same procedure as that used above for Pa, 
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namely by the mass renormalization counterterm. The result obtained at the end of this 
procedure is again proportional to m and vanishes when m —>■ 0. 

Note that the combination n(0,0;m, T) — no(0,0;Mo) in Eq. ( |175| ) enters also the cal- 
culation of the renormalized quark-quark scattering amplitude to order g"^ in perturbation 
theory: 

T(0, 0) = / - [n(0, 0; T) - no(0, 0; Mo)] = / - 7 In -^^^^^ (176) 



where the limit m — > has already been taken (we used Eqs. (|A2|) , ([A4|) and ( |A18|) ). It can 



easily be verified that this expression is invariant under renormalization group transforma- 
tions (to order g^) when g runs according to the lowest order /^-function (0). Now, as usual 
in perturbation theory, one can choose the scale Mq so as to minimize higher order correc- 
tions. In the present case, the order g^ contribution to the scattering amplitude vanishes for 
Mo = 7rT/e'^-E+i_ Note that for this choice of Mo, 

/(M„ = .r/e™«) = ^. (177) 

which, is the coupling constant g^ introduced above, in Eq. ( |154| ). Note also that, for this 
choice, the quark-quark scattering amplitude in Eq. ( |176| ) coincides with that calculated 
within the non-perturbative scheme (see Eq. (|149|) ). 

Before we move on to the last order g^ contribution (the only one to survive the limit 
m —>■ which we implicitly assume from now on), let us return to the contribution Pf, given 
in Eq. ( |170|) , in order to make closer contact with the calculation. One can rewrite Pb 
as 

-ioo+e 

(By performing the w-integral and suitably rearranging the statistical factors, one can re- 
cover, form this expression, Eq. ( |171| ) above.) The quantity Il{uj,q]T) is ill defined, since it 
contains an infinite constant contribution, independent of the temperature. We can absorb 
this constant into no(0,0;Mo) and discard the corresponding vacuum contribution. Next 
we write: 

n(cu, q; T) - no(0, 0; Mq) = [n(u;, g; T) - n(0, 0; T)] + [n(0, 0; T) - no(0, 0; Mo)] . (179) 

The last term within brackets is (1/vr) In (7rT/Moe'''^+-^) (see Eq. ( |1 76| ) above). The first term 
within brackets is {1/2'k)Bt{Q/T) which, at large Q, behaves as {l/2T^)\n{Q'^T'^ / {T^M})) 
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Pb = ^ r''^ ^ (1 + 2N{u;)) Uiu, q- T) = n(r = 0, a; = 0; T). (178) 



(see App. A). This suggests the following temperature independent subtraction which trans- 
forms Ph into 



2 ^ 
2^ 



(180) 



At this stage, it is necessary to separate the integral into the part which contains the 
statistical factor, and the other which does not. The part with the statistical factor is 



calculated by deforming the contour in the usual way (see subsection [IV Bj) and yields 



/j2 /.oo POO 

Ph=^^l I dujN{uj)lmBT{oo/T,q/T). (181) 

This contribution is independent of the vacuum subtraction. Changing to the dimensionless 
variables s = uj/T, r = q/T, fixing the scale Mq at the value Mq = ttT /e""^^^ and identifying 
(yf^ with the running coupling at that scale (i.e., replacing by g'^ in Eq. ( |181| )), and finally 
multiplying by (47r/M|) (see Eq. (|T26|)) and dividing by {T^ /M]){g^ /27i) (see Eq. ([T62|)), 
one recovers the expression ( |163| ) of 6*2- The part without the statistical factor in PH reads 

and it can be verified that it is finite. By choosing again Mq = vrT/e'^^"'"^, and substituting 
g"^ for g^ in Eq. ( [L82| ), one gets: 

The same manipulations as done before for P^ allow us to recover the expression (|156|) of 

F2. 

At this point, we have all the necessary ingredients to understand the origin of the 
cancellation F2 + G2 = observed in the previous subsection. We have just seen that 
Pb = P^ + Pff oc F2 + G2; but from the calculation of Pb done above with Eq. ( |17(JD we 
have Pb = 0. Hence the cancellation. Consider finally 

P, = ^J[d'Q][U{u;,q;T)f. (184) 

We can proceed as we did for Pc and first combine Pd with the contribution coming from 
the coupling constant renormalization of Pb which transforms Pb into P^': 

Pi = Pb{l + g^UoiO, 0; Mq)) = Pb - ^ J [d^QMu, q; T)Uo{0, 0; Mq). (185) 
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(Of course, according to the previous discussion, this is a vanishing correction, but it is 
convenient to carry it along as it makes the procedure more transparent.) By adding to Pd 
the order g'^ contribution from P^', we get: 



P' 



^ I [d^Q] [Uico, q- T) - no(0, 0; M,)f 



(186) 



where we have also added an irrelevant constant ~ (no(0, 0; Mq))^. The divergence in 



Eq. (|186|) is dealt with in analogy with what we just did for P^- Thus we write 



d^g 
4 7 (27r)2 



Q 



27r^''(,M2e2 



(187) 



and it can be verified that P^' is finite. 



At this stage, one may choose the scale to be Mq = ttT/q"^^'^. Then P^' simplifies into: 



~9' 



d2g 



27r 



(188) 



4 J (27r)2 

We can now perform the cu-integral as we did earlier and separate the contributions contain- 
ing one or no statistical factors. We obtain in this way the contribution of order to the 
functions F and G. These are: 
1 



P4 = — 



TT JO 



Ga 



da dr I [Briis, r)f - [in (r^VAf|(s^ + 



-1 Td. 



drA^(s)ImPr(s, r)ReP5n(s, r) 



(189) 
(190) 



which agrees with the expressions given above (Eqs. ( p.57| ) and ( |164|) ). 

To summarize, the perturbative calculation performed in this last subsection has allowed 
us to recover the high temperature limit obtained previously within the 1/A^ expansion. 
The calculation involves a running coupling at a scale Mq which can be chosen, as usual in 
perturbation theory, so as to minimize the high order corrections. Here the choice Mq = 
T^T/e""^'^^ allows us to identify the running coupling with g which naturally emerges in the 
high temperature expansion of the order 1/A^ calculation. We have seen that the second 
order terms vanish, because they are proportional to the quark condensate which vanish in 
the high temperature phase; this explains the cancellation of F2 + G2 shown in the previous 
subsection. Finally, we have also been able to reproduce explicitly the analytical formulae 
giving the first few terms of the high temperature expansion of the pressure. This analysis 
confirms that, as could be expected from asymptotic freedom, the behavior of the system 
at high temperature can be understood in terms of perturbation theory. 
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VI. CONCLUSIONS 



The Gross- Neveu model in 1+1 dimension has proven to be a very useful playground to 
study various aspects of renormalization at finite temperature. We have obtained, at next- 
to- leading order in the 1/A^ expansion, a new expression for the effective potential, which 
is explicitly invariant under renormalization group transformations, both at zero and at 
finite temperature. We have verified, in a non-perturbative context, that the temperature 
dependent ultraviolet divergences cancel, as expected from general arguments: the same 
counterterms which make finite the effective potential at zero temperature also make it finite 
at non zero temperature. There is presently a lot of interest in exploring non perturbative 
techniques to study the thermodynamics of gauge fields, in particular of QCD ||3^, ^ ^ 
and the present investigations contribute to this general effort, in a much simpler context 
than that of gauge theories. In contrast to other methods, such as for instance those based 
on $ derivable approximations |]39| , the 1/A^ expansion has the advantage of being at the 
same time non perturbative and of providing an expansion parameter: the latter is helpful 
to understand the systematics of the cancellation of ultraviolet divergences. Within this 
framework, there are obvious generalizations of the model studied in this paper, for example 
considering continuous symmetries and higher dimensions, which would be worth exploring. 

Part of the motivation for studying the Gross-Neveu model specifically in 1+1 dimension 
is that it shares with three dimensional non abelian gauge theories like QCD the property of 
asymptotic freedom. This leads us to expect that at high temperature, the system behaves 
as a weakly interacting gas of the original constituents, somewhat analogous to the quark- 
gluon plasma of QCD. The detailed analysis of the high temperature limit is an important 
part of the present study. We have verified that, as expected from asymptotic freedom, the 
pressure can be expanded in terms of an effective coupling which decreases with increasing 
temperature. In fact, we have shown explicitly that, in this high temperature regime the 
prediction of the 1/A^ expansion is identical to that of ordinary perturbation theory. We 
find that the pressure goes slowly towards that of the ideal gas as the temperature increases. 
Compared to QCD where a similar phenomenon occurs, one finds that in the present calcu- 
lation the pressure is higher than that of the free gas, while it is lower in QCD. In QCD, the 
dominant correction to the free gas behavior is of second order in the coupling strength, and 
the sign of the effect is easily understood from the corresponding expression of the entropy 
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density ||4^. Here the contribution of the order g"^ cancels and the leading correction to the 
free gas behavior is of order g'^. Even if the basic mechanisms at work are different here and 
in QCD, it would be interesting to explore further the properties of the high temperature 
phase of the Gross-Neveu model along the lines developed in Ref. P0|. 

We find that the expansion provides a sensible approximation scheme at both low 
and high temperature: At low temperature, this is because it produces only small correc- 
tions to the mean field physics; at high temperature, it reproduces essentially perturbative 
calculations. But, for temperatures of the order of the mean field transition temperature 
the approximation breaks down: a Landau pole appears in the calculation and fluctuations 
become too large to be treated as corrections. (The Landau pole is harmless at low temper- 
ature and disappears at high temperature.) It would be interesting to explore further this 
breakdown of the expansion, in particular to study to which extent it is a consequence 
of the one-dimensional character of the system, by using different techniques, such as Exact 
Renormalization Group Equations, with which fermionic models have already been exten- 
sively studied at zero and finite temperature |^ (for a review see [0). In this context, 
bosonic fluctuations can be included in a simple way, but the 1/N expansion has not been 
explored beyond the leading term. 
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APPENDIX A: RENORMALIZED a PROPAGATOR 



The bare fermion loop is defined in Eq. (p2D and reads: 

+ K -K' 



n(Q;M,T) = 2 j{d^K}- 



(Al) 



'~K'^ + M^){-K'^ + M^) 

where the notation / {dpK} is that of Eq. and K' = K + Q. In this equation, and 
throughout this Appendix, M is considered as a given parameter, not necessarily equal to 
the physical fermion mass Mf. Note that Il{Q) is dimensionless. 

One may perform the sum over the Matsubara frequencies (see Eq. ( |X9| ) below) to verify 
that n can be written as a sum of a zero temperature contribution Ho and a finite tem- 
perature contribution 11 which goes to zero when T 0. This separation is well defined 
only when M ^ 0: indeed the zero temperature contribution, IIq{Q,M) is divergent when 
M — > 0. At finite temperature, the fermionic Matsubara frequencies provide an infrared 
cut-off which makes Il{Q; M 0,T) well defined. 

We consider first the zero temperature contribution IIq{Q] M). This may be obtained 
directly from Eq. ( [Al| ) as an Euclidean integral, function of Q|; = Qq + only. Once 
regularized with a cut-off A, it yields: 



Iio{QE-:M) = — 



(A2) 



where 



Bo(z) 



14: + Z 



In 



V4 + 2; + 



(A3) 



+ Z - y/z 

is an analytic function of z with a cut on the negative real axis. Some properties of Bo{z 
for z = X real and positive are useful. For small x 

n2 



while for large x 



Bo{x) 



B,{x)=2 + l-^- + 0{x^ 



1 2 f a + b\nx 
Inx + -(1 + Inx) + O 



(A4) 



(A5) 



The expression 1 + q'^UqIQ; M) may be interpreted as the bare inverse propagator for 
the a field in leading order of the expansion (see Eq. (pT]) ). The renormalized inverse 
propagator is obtained by multiplying it by Z (coming from the renormalization of a in 
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Eq. (HI)), and replacing gl by g^/Z, with Z = = (g^/Ti) ln(A/M/) (see Eq. (|TD with 
Z^^^ = 0). One then obtains: 



Do\Qe;M) 



z(l + '\,iQ,;M))^f^ 



In 



Bn 



M2 



9\ M\ g' 1 ^ 



Q 



4M2 



1 + W1 + 



4M2 



(A6) 



This is an increasing function of Ql;, whose minimum, at Ql; = 0, coincides, up to a factor g"^ 
(see Eq. (p7|)), with the second derivative with respect to M of the leading order contribution 
to the renormalized effective potential (Eq. (|55D): 



0; M = — 1 + In 



dM2 ■ 



(AT) 



Note that D^\Qe = 0; M) vanishes at the value M, = M//e. For M > M^, Dq^{Qe; M) 
is always positive; for M < M^, Dq^{Qe; M) vanishes at some positive value of Ql; and 
is negative for smaller values of Q|;. That is, when M < the a propagator has an 
unphysical pole, whose role in the next-to-leading order calculation of the effective potential 
is discussed in the main text. 

For large values of Q|;, we obtain from Eq. 



2tx 



In 



+ 



2M2 



Q 



1 + ln 



M2 



+ 



M2 



(A8) 



We now turn to the finite temperature contribution, n(cc;, g; M, T). This is obtained from 
Eq. ( |A1| ) by performing the sum over Matsubara frequencies. The resulting expression for 
the total n is: 



no(cu,g;M) + n(cu,g; M,T) = 
^_ + M2 - fc2 + g74 



dfc 

-oo 27r 
n^ — n_ 



n< — 



+ 



\uo- E+ + E_ uj + E+-E_J 
q'^/4: / + n_ — 1 n_|_ -f- n_ — 1 



2E,E_ 



io-E, 



E_ uj + E+ + E_ 



where 



E. 



± 



{k ± g/2)2 + M2, 



n± = n{E±) 



1 



(A9) 



(AlO) 



e/3i?± + 1 • 

The finite temperature contribution proper, i.e., n(to', g; M, T), is obtained by keeping 
only terms proportional to statistical factors (i.e., by dropping the I's in the third line 
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of Eq. (|A9|) ) . One can then verify that 



fl{uj, q) = n(-cj, q) = n(u;, -q). 



(All) 



Thus we need to calculate 11 only in the region uj > 0, q > 0. Note also that changing k 
into —k in the integrand of Eq. ( [A 91 ) is equivalent to changing q into —q which leaves the 
integrand invariant; we can therefore limit the /c-integration to the range [0, oo[. 

The imaginary part is defined as usual as the imaginary part of n(ti; + ie, q) with uj real. 
In the region of interest, this is non vanishing for k values such that (a): — = u, or 
(b): E+ + E_ = UJ. The case (a) is typical of finite temperature and corresponds to scattering 
processes in the heat bath. The case (b) corresponds to processes which are analogous to 
those occurring at zero temperature (decay of the sigma into quark-antiquark pairs; see 
Sect. |VB|) . For each value of uj and g, the "phase-space" is limited to a single value of the 



momentum and the energy of the quark. This is given by {ujq = \J(f ^ AM'^): 



< UJ < q 



(scattering processes) 



uj'^ — uj"^ 

— Uj"^ 



El 



uj"^ — Uj"^ 
q^ — uj"^ 



(A12) 



UJ > UJ„ 



(pair creation) 



UJ^ 



UJt 



k" 



E' 



2 ^ ^2 - g2 

The corresponding imaginary parts are given by: 



- \ UJ ± q, 
2 ' ^ 



(a) 



<uj < q 



imn 



2\ q 



UJq — Uj"^ 



UJ^ 



UJ^ 



UJ:, 



\ Cj2 - g2 



n{El)-n{El] 



(A13) 



(A14) 



UJ > UJa 



Imn = 



ct;2 — itjl 



n{EX) + n{E^_ 



(A15) 



2^ w2-g2 

For u; > outside the regions [0, g] and [tu^, cxd[, Imll = 0. The real parts are determined as 
principal value integrals. 

Note that the contribution of the scattering processes to the imaginary part is negative, 
while that of the pair creation processes is positive. In fact, in order to get a physically 
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meaningfull result, the latter should be added to the (negative) zero temperature contri- 
bution (see Eq. ( |142| ) ) : the finite temperature contribution can then be interpreted as a 
"suppression" of pair creation due to Pauli blocking. 

Consider now IT in Eq. ( |A9| ) , and take the limit M —>■ 0. We get: 

(A16) 



n(a;,g; T) = 
where we have used the property 



A/2 dk 4k 

TT CU^ — 



1 



1 



3/3(fc+g/2) _^ I Q(3{q/2-k) ^ ]_ 



1 



(A17) 



e/3(fc-g/2) + 1 6^(9/2-^) + 1 

to rearrange terms. The ultraviolet cut-off A/2 (rather than A) is chosen so that Z + g^U 
coincides with the previously defined renormalized inverse propagator in the limit M — > 0. 
Indeed, we have 



, /•A/2 dk 1 



1 



1 



(eT^^A) 



TT 



(A18) 



where the second equality holds up to terms which vanish when A ^ oo. Combining this 
result with the expression ([51| ) of the renormalization constant Z (with Z^^^ = 0), and using 
Eq. (|n|), we get 

.2 



Z + g'U{0,0;T) 



9 



In , 



(A19) 



As expected, up to the factor g"^, this is the second derivative of the effective potential at 
M = (see Eqs. and (p9|) ). 

The inverse propagator at M = is then: 

D-\uj, q;T) = Z + g'U{uj, q;T) = Z + g'U{0, 0; T) + g'[U{uj, q; T) - n(0, 0; T)] (A20) 

where the combination n(ci;, q; T) — 11(0, 0; T) = {1/2'k)Bt{Q/T) is finite. Thus, in analogy 
with Eq. ( |141|) , we write 



D-\uj,q-T) = ^ 



277 



(A21) 



Thus defined, the dimensionless function Bfi^J /T^q/T) behaves as ln(T^^/M|(— cj^ + 
g2)/T2) + C(T^/(g2 _ ^j'^f) when g > T or a; > T. At small a;,g, 5t is regular 
{Bt{u}/T, q/T) ^ as cu, g 0). 
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FIG. 17: The function Re5r(s,r) 

The imaginary part of D^^iuj, q; T) is finite and given simply by {u real and positive): 

lmU(LJ,q]T) = - \- \- V (A22) 

At small uj, and as long as g ^ T, Imn(ct;, g; T) is a linear function of uj, with a slope 
—q/AT. When q > T however, the slope vanishes exponentially, as (q/T) exp{—q/T). For 
large values of q/T, the imaginary part becomes a step function of height -1/2 located at 
u = q. 

The real part is obtained numerically from the principal value of the integral (|A16 ), 
by combining with the renormalization constant to eliminate the cut-off. It can be also 
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expressed in terms of the real part of the function Bt'- 

X 



ReBris, r) = PP / dx 
Jo 



x^ — X 



(A23) 



where we have set x = k/T, s = uj/T, r = q/T, and 

1 



/: 



± 



+ 1 



e2 



(A24) 



A plot of ReBT^s, r) is given in Fig. 0. 

Let us finally examine the asymptotic behavior of Il{uj, q) when q (or uj) is large. This is 
dominated by the zero temperature contribution and the leading term is given by Eq. ( [A 81) . 
However, at next-to-leading order, there is a relevant contribution coming from 11. To 
estimate it, we go back to Eq. (|A9|) , change u — > iqQ and to go to polar coordinates: 
go = rcosO, q = rsin6'. Next, one performs changes of variables such that all the statistical 
factors in Eq. (^3) have the same energy argument. Then it is not hard to get: 



mo, q)- I [8^^ - 8 sin^ OiEl + P)] 1 



In particular. 



/•2^d^~, 
Jo 2n 



4M2 



=° dk Uk 
oo 2tt Ek 



(A25) 



(A26) 



APPENDIX B: THE FUNCTIONS Fo(x), F{x,T) AND G{x,T) 



The function Fo{x) introduced in Eq. ( [73D of the main text is defined by: 



I roo ( 



Inx + Bq{u/x) 



X — 1 X — 1 — xlnx 

2 2 

u u\nu 

Inx + Bq{u/x) 



2 Jo 



du In 



B.iu) 



(Bl) 



The last two terms in the first line of this equation correspond to the A-dependent terms in 
Eq. (|65D, which we have rewritten as: 

1 A2 _ /-A I d^Q 1 
A^^Jd] ~ JeMf (27r)2 Q2 + ^' 
/•A 1 rl'^n 

(B2) 



In In 



47r 



M 



f 



eM 



, (27r)2g2in(g2/M2)- 
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Thus, by combining the integrand as we did in Eq. (|B1|), we have obtained convergent 
integrals and the hmit A ^ oo could be taken. Note that, as stated in subsection |111 B| , the 
function Fq{x) is defined only for x > Xm, = 

To find the first derivative Fq, it is convenient to separate the contribution of the last two 
terms of the first line of Eq. (|B1|), which we call -Fct- Since F^rj, is divergent, we re-introduce 
a cut-off A, and get: 



In order to calculate the remaining contribution to Fq{x), we use the change of variables 
to get 

d?; In 11 

F;,{x) - F^^ix) = / ^Zi, ^ (B5) 

Ji V inx + In V 

v—l 

where a = v{u = A"^) ~ A^/x for large A. In order to recombine both contributions (B3) 



and ( |B5|) , and let A — ^ oo, one observes that 



— -1 + - =/ — + - — -lnx + C(lnx/lnA2). (B6) 

2 u y mu J Je2 u \ mu J 

Thus, adding ( P5| ) and ( PB| ) and now letting A — > oo, one gets after simple algebra: 

du / 1 (u- 1)^ \ 



F'ix) =\nx T — 

Je2 U 



«2 



e2 u \\nu (m^ — 1) Inx + (m + 1)^ Inw^ 
du In u 



U— 1 

The second derivative Fq follows easily from the previous expression. One has: 

du Inu 



r°° du I I Inu \ f 

° + 1 = ^ ~ []^ ~ (inx+^ln«)2 )~ Ji 



1 u (In X + In u 



(B8) 



Finally, the third derivative Fq", needed for instance to calculate the /5-function at next-to- 
leading order, is obtained by taking the derivative of the expression above: 



oo 



du In u 



x'F'"(x) + xF''(x) = 2 / — — y- (B9) 
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Note that F^, xFq and x'^Fq + xFq are actually functions of the variable Inx. In the main 
text we use the limit of these expressions as In a; ^ oo. This limit is easily obtained from 
the equations above. On finds: 



Fl^{x) ~ (lnx)(lnlna;) xF^{x) ~ Inlnx and xF^ + x^Fll'{x) ~ 1/lnx. (BIO) 

We turn now to the function -F(x, T) which is obtained when we regroup what remains 
of V^^i after elimination of the temperature dependent divergence in Eq. (123) with the 
function Fq introduced above, Eq. (pl|) . We define: 

'D-^{r,e]M) 



-Lf{xm,T) = -^Fo{xm) 
An Air 



1 r'^f rdr /•27r 



1 rdr I d9 ^ 



2 Jo 27r 7o 27r 
L)-i(r,^;M) 



In 



DQ\r^;M 



f) 



+ 



8M^ 



du- 



r2 1n(r2/M2) \Jo ^u^ + Xm 



(Bll) 



In this equation we have set lo = iqo and changed to polar coordinates: go = rcos6, 



q = rsin^. Also, we use abusively the notation D ^{r,6) for D 



UJ 



ir cos 6,q = r sin i 



and similarly for Dq^. To get the last term in the expression above we used the second 
equation in (|B2D and the fact that InlnA/M/ = InlnA/M + C(l/ ln(A/M/)). Note finally 
that, since F contributes to a term of order we have set xm = M'^/Mj, ignoring the 
factor of order in Eq. (P^). One has: 

r^/^de 



F{x,T) 



2v dv 







TT 



\n{A{v,e;x,T)) 



X — 1 



+ 



2v dv 
Jo 



f 2 In f 2 

/2 de 

IT 



l^a; — 1 — X In X — Ax 
\n{A{v,e;x,T)) 



du—z. 

yu^ + X . 



(B12) 



where the function A{v, 6; x, T) is 



A 



In X + Bo{v^/x) + 27ifl{v, 9; x, T) 



(B13) 



In ( P12|) we have used the symmetry properties ( |A11|) to limit the integration over 6 to the 
interval [0,7r/2]. Note finally that F{x,T = 0) = Fo{x), with Fq given by Eq. (0). The 
function F{x,T), evaluated numerically, is shown in Fig. |18|, for various values of T. Note 
that, when T < Tc, F{xm,T) is defined only for xm > xmS"^) (see subsection |IV B[). 



70 



-3 



FIG. 18: The function F{x,T) for various temperatures. Pull line: T = 0, i.e., Fo{x); dotted line: 
T = 0.565M/ = Tc, dashed line: T = Mf 

To study the high temperature limit of the the function F{xm,T) (needed in subsection 
VD|) it is useful to write F{xm = 0, T) in a different form. Adding and subtracting the 



expression 



Jo Jo 71 \\n(VTT^) ^ ' 



we can rewrite F{0,T) as 

ln(TVT2) + r) 



+ (s2 + r2)TVM 



oo 



+ / 2wrft; W — In ^— + ^ + ^ > + / 2t;rft; / — In 



IT \ Bq{v'^) J t>2 ti^lnf^j Jo JO TT Y Bq{v'^) J 

(B15) 

In the first line we introduced the function Bt (defined in App. A) and made the change of 
variables s = {Mf /T)v cos6, r = [Mf /T)v sinO. It is a finite expression (see the comment 
after Eq. (|A21|) ). The second line is also finite (see Eq. (|A5|) ) and it is independent of the 



temperature. 

Finally the function G is 



A POO rOO 

G{x,T) = / dr dsN{s)S{s,r-x,T). (B16) 
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FIG. 19: The function G(x, T) versus x for various temperatures T (in units of the fermion mass 
Mf). The dotted Une corresponds to T /Mf = 0.1 and indicates the typical behavior of the function 
for smaU values of x when T < Tc- For the other temperatures the function is plotted only for 

X > XM-AT)- 

where 6 is given by Eq. ( |129| ). We have set here s = u/Mf, r = q/Mf and N{s) = l/(e^'^ — 1). 

The function G{x, T), evaluated numerically, is presented in Fig. [1^ for various temperatures. 

One can understand the behavior of G at small xm from the behavior of the phase 

shift. At moderate temperatures it is enough to consider the zero temperature phase-shift 

in evaluating the integral in Eq. ( |B16| ). For M = Mf, using the expression of the zero 

1-1/ 



(B17) 



temperature propagator Dq \u, q; Mf) given in App. A, we find {ujg = \Jq^ + 4Mj): 



tan(5) 



TT 



2 Arctanhy^ (u;^ — u'^) / (cj^ — q^) 
Consider this expression for fixed g, as a function of u. For u ojq this behaves as q) 



7r/2 - (2/7r)^(c<;2-cu2)/(cj2-g2) ^hile at large 6{uj,q) ~ 7r/\n{uj^ - q^). For M ^ My, 
DQ^{uj,q]M) = DQ^{u,q;Mf) + (^Vtt) ln(M/M/). A plot of ReDQ^LU , q; M) is given in 
Fig. 1^. For M < Mf, the real part of DQ^(uj,q;M) vanishes at points cUa and Ub and is 
minimum at the point Ug where it vanishes for M = Mf {ua < ojq < oob). For u < Ua, 5 is 
small and positive; it jumps to a value close to vr when uj crosses the value uja-, and remains 
approximately constant until ~ cjg, at which point it decreases rapidly to reach the value 
— 7r/2 at ujb] it then goes slowly to zero as a; ^ oo. When M M^{T) from above, Ua 0, 
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and the region uj < Uq, where 5 ~ vr, gives a large contribution to G, amphfied by the large 
value of the statistical factor at small uj, where A^(ci;) ~ T/uj. This is the origin of the large 
negative contribution at small xm that is visible in Fig. [1^, particularly in the curve labelled 
T = O.lMf. 



APPENDIX C: THE FERMION SELF-ENERGY 

Here we calculate the fermion self-energy defined in Eq. (|35D . We have: 

E(A-) = -f^ / [d'Ql M) JfJ^Zg/Zp " + 

where M is the renormalized mass, and g the renormalized coupling. The functions a{K) 
and b{K) can be obtained easily by projection: 

= ^.-^^mK) = --4 /id^Ol D(Q; M) (C2) 

1 n'^ r M 

KK) = ) = jl^^Q] DiQ; M) _ ^ (C3) 

One can separate the self-energy S into a zero temperature contribution and a finite 
temperature one, as we did for the fermion loop in App. A. Accordingly we write: 

= So(i^) + t{K) (C4) 

and similarly for the functions a{K) and h{k). It will be verified shortly that a{K) is finite, 
while h{K) is divergent. The divergence is a mass correction which is eliminated by the mass 
counterterm Introducing a cut-off A, one can write: 

M /AM M , 
KA-) = --lnl„(^]+-iW. (C5) 

where, as we shall verify, h'{K) is finite. The full fermion two-point function takes then the 
form: 



S-\K) = -^(l-a) + M(^l + ^j 



b{K) 



= -${!- a{K)) +m(i + ^ + b'{K)^ (C6) 

where we have used Mb = TZ'M, = 1 + Z''^^^ /2N and Z'^^) = lnln(A/M/) + Z'^'^\ 
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We now proceed to the explicit calculation of the zero temperature contribution to S. 
The corresponding functions ao and Bq are functions of only. By going over to Euclidean 
momenta, and performing the angular integrals, one obtains: 



b'oiKj 



9 



A2 



dQ 



DoiQhM) 



(C7) 



and 



3_ 

'n 



Do{Ql;M)\l + 



One can verify on these explicit expressions that both ao(-ft'^) and Bq^K'^) are finite. 

In fact we shall need these functions only for M = Mf and Kj^ = 
ao = aoiKl = -M|), b'^ = h'^{Kl = -M]). We have: 



(a 



M|, and we shall set 



ao 



Do{Ql;Mf 



N Jo 2n (Q2 + 2M])JQ^ + 4MfQ' + + AQ'^M\ 



dQ' 



(C9) 



(CIO) 



2N Je^Mj g2 ln(Q2/M2) N Jo ill ^Q^ + 4M]Q^' 
By using the explicit expression ([A6|) of the propagator together with the change of variables: 



u 



Q' 



+ 



Q' 



+ 1, 



\ iMj \ AM] 

one can simplifies these expressions, and get: 

dw — 1 1 



1 u"^ + 1 Inw^ 



, du - 1 1 r"^ du 2 1 
b'^ = - I r ; + 



(Cll) 



(C12) 



(C13) 



'i M + 1 Inu^ Je u u' + 1 Inu"^ 
We are now in a position to calculate the change in the fermion mass coming from S. 
We have, in the vicinity of the mass shell, to within terms of order 1/N': 



= -M/(l - a) + Mra^ U + ^ + ^' J 



(C14) 



where we have replaced M by the value Mmin it takes at the minimum of the effective 



potential. Writing = M^m + and ignoring terms of order we get: 

a + h' + — — = -— + ^ 



M 



2N 2N N' 
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(C15) 



where 



°° du - 1 1 




du - 1)2 1 



(C16) 



e it^ rt^ + 1 In u' 



u'^ + 1 In 



Numerical integration gives (p = —0.126229. 

Finally, we use in the main text the expression of the propagator in coordinate space. 
This is easily obtained from the following "mixed" representation 



where A± = (Ep ± 70 (7V + ^)) /S-Bp and Ep — ^/p'^ + m'^. 



[1] D. J. Gross and A. Neveu, Phys. Rev. D 10, 3235 (1974). 

[2] B. Rosenstein, B. J. Warr and S. H. Park, Phys. Rev. Lett. 62, 1433 (1989). 

[3] G. Parisi, Nucl. Phys. BlOO, 368 (1975). 

[4] K. Shizuya, Phys. Rev. d21, 2327 (1980). 

[5] J. Zinn- Justin, Nucl. Phys. B 367, 105 (1991). 

[6] D. J. Gross, In Les Houches 1975, Proceedings, Methods In Field Theory, (Amsterdam, 1976, 
141-250). 

[7] J.C. Conins, Renormalization, (Cambridge University Press, 1984). 

[8] N. P. Landsman and C. G. van Weert, Phys. Rept. 145, 141 (1987). 

[9] M. Le Behac, Thermal field theory, Cambridge Monographs in Mathematical Physics, (Cam- 
bridge University Press, 1996). 

[10] R. F. Dashen, S. k. Ma and R. Rajaraman, Phys. Rev. D 11, 1499 (1975). 

[11] L.D. Landau and E.M. Lifshitz, Statistical Physics (Pergamon, London, England, 1958). 

[12] P. Fendley and H. Saleur, Phys. Rev. D 65, 025001 (2002) [arXiv:hep-th/0105148]. 

[13] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 124, 246 (1961). 

[14] L. Jacobs, Phys. Rev. D 10, 3956 (1974). 

[15] B. J. Harrington and A. Yildiz, Phys. Rev. D 11, 779 (1975). 

[16] W. Dittrich and B. G. Englert, Nucl. Phys. B 179, 85 (1981). 

[17] U. Wolff, Phys. Lett. B 157, 303 (1985). 



S{t >Q,p) = A+7oe-^-"(l - Up) + A_7oe^-"n, 



'p 



S{t<0,p) 



A+7oe-^^^np - A_7oe^''^(l - Up) 



(C17) 



75 



[18] A. Barducci, R. Casalbuoni, M. Modugno, G. Pettini and R. Gatto, Phys. Rev. D 51, 3042 
(1995) [arXiv:hep-th/9406117]. 
3] V. Bernard, U. G. Meissner and I. Zahed, Phys. Rev. D 36, 819 (1987). 
D] T. Hatsuda and T. Kunihiro, Phys. Lett. B 185, 304 (1987). 
1] M. Asakawa and K. Yazaki, Nucl. Phys. A 504, 668 (1989). 

2] A. Barducci, R. Casalbuoni, M. Modugno, G. Pettini and R. Gatto, Phys. Rev. D 55, 2247 

(1997) [arXiv:hep-ph/9607457]. 
3] M. Modugno, Riv. Nuovo Cim. 23N5, 1 (2000). 

4] J. Hufner, S. P. Klevansky, P. Zhuang and H. Voss, Annals Phys. 234, 225 (1994). 
5] W. Florkowski and W. Broniowski, Phys. Lett. B 386, 62 (1996) [arXiv:hep-ph/9605315]. 
3] R. G. Root, Phys. Rev. D 10, 3322 (1974). 
7] R. G. Root, Phys. Rev. D 11, 831 (1975). 
3] J. F. Schonfeld, Nucl. Phys. B 95, 148 (1975). 
3] R. W. Haymaker and F. Cooper, Phys. Rev. D 19, 562 (1979). 

D] J. P. Blaizot and G. Ripka, Quantum Theory of Finite Systems, (MIT Press, Cambridge, 1986). 
1] J.L Kapusta, Finite temperature field theory, Cambridge Monographs in Mathematical 

Physics, (Cambridge University Press, 1989). 
2] L. S. Brown, Quantum, Field Theory, (Cambridge, UK: Univ. Pr., 1992, 542 p). 
3] J. Zinn-Justin, Quantum Field Theory And Critical Phenomena, International series of mono- 
graphs on physics, 77, (Oxford, UK: Clarendon, 1989, 914 p). 
4] G. Ripka, Quarks Bound By Chiral Fields: The Quark- Structure Of The Vacuum And Of 

Light Mesons And Baryons, (Oxford, UK: Clarendon Pr., 1997, 205 p). 
5] S. R. Coleman and E. Weinberg, Phys. Rev. D 7, 1888 (1973). 
3] W. E. Thirring, Annals Phys. 3, 91 (1958). 
7] W. Wetzel, Phys. Lett. B 153, 297 (1985). 
B] L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 

3] H. van Hees and J. Knoll, Phys. Rev. D 65, 025010 (2002); H. Van Hees and J. Knoll, 

arXiv:hep-ph/0111193. 

D] J. P. Blaizot, E. lancu and A. Rebhan, Phys. Rev. D 63, 065003 (2001) [arXiv:hep- 
ph/0005003]. 

1] J. O. Andersen, E. Braaten and M. Strickland, Phys. Rev. D 63, 105008 (2001) [arXiv:hep- 



76 



ph/0007159]. 

[42] U. Ellwanger and C. Wetterich, Nucl. Phys. B423, 137 (1994). 

[43] J. Bergcs, D.-U. Jungnickel and C. Wetterich, Phys. Rev. D59, 034010 (1999). 

[44] J. Berges, N. Tetradis and C. Wetterich, Phys. Kept. 363, 223-386 (2002). 



77 



